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PROGRESS REPORT 

This report summarizes research carried out at Harvard 
University under NASA grant NGL 22-007-247, during the report 
periods 1 August 1975'-31 January 1975 and 1 B’ehruary 1975^31 July 
1975- Pour papers and abstracts have been published since the 
last such report (1 August 197^) and five other papers have been 
accepted for publication and are now in press. These publications 
are listed here and copies of some of them are appended. 

Research during the report year has followed a number of 
different lines. Our cooling-rate studies of 12002 have now been 
completed and our results have been interpreted In terms of the 
crystallization history of this rock and certain other plcrltlc 
Apollo 12 samples. This work was reported at the Sixth Lunar 
Science Conference and a paper on the experimental portion of the 
study has been accepted for publication in the Geological Society 
of America Bulletin. Calculations of liquid densities and vis- 
cosities during crystallization, crystal settling velocities, 
and heat-loss by the parent rock body are in progress, as are 
petrographic studies of other Apollo 12 samples with the aim of 
generalizing and extending the results of the experimental study. 

We have given some thought to the process of magmatic differ- 
entiation that must have accompanied the early melting and chemical 
fractionation of the moon ^s outer layers. The effects of the 
pressure differential in such an immense magma body Introduce some 



novel features into the crystallization process that seem not to 
have been previously considered. The chemical and petrological 
nature of the outer few hundred kilometers of the moon and its 
potential as a source region for mare basalts cannot properly be 
understood without taking these features into account. This work 
was also reported at the Sixth Lunar Science Conference and a 
paper has been accepted for publication in the Proceedings volumes 

Our Interest in the source regions of both high- and low- 
titanium mare basalts continues. The work on hlgh-Tl basalts 
reported last year has been widely cited and used as a basis for 
geochemical calculations by a number of authors, but it has also 
been criticized and challenged, particularly by the group at 
Canberra. We have had many discussions on this problem and 
expect to conduct further experiments on a critical saimple, 7^275, 
in an attempt to resolve the points of controversy. We are also 
preparing a manuscript on the nature and origin of the low- 
titanium mare basalts. 

Finally, we are attempting to analyze a very large body of 
chemical data on our experimental run products and to summarize 
this data in forms useful to terrestrial petrologlsts and geo- 
chemists as well as to other lunar scientists. A first report 
on this work was given at the annual meeting of the American 
Geophysical Union. 



PUBLICATIONS 

(1 August 197^ - 31 July 1975) 

The petrology of the Apollo 17 mare basalts. J. Longhl, D. 
Walker, T. L. Grove, E. M, Stolper, and J. F. Hays. Proc. 
Fifth Lunar Sc i. Conf., Geochlm. Cosmochlm. Acta, Suppl. 5> 

I, 447-469, Pergamon Press, 1974. 

Differentiation of a very thick magma body. D. Walker, J. 
Longhl, and J. F. Hays. Lunar Science VI , 838-840, The Lunar 
Science Institute, 1975. 

12002 Revisited. D. Walker, R. J. Kirkpatrick, J. Longhi, 
and J. F. Hays. Lunar Science VI , 841-843, The Lunar Science 
Institute, 1975. 

Pe-Mg distribution between olivine and lunar basaltic liquids 

J. Longhi, D. Walker, and J. F, Hays. Trans. Am. Geophys. 

Un. (E®S) 56, 471, 1975. 
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PAPERS IN PRESS 

1. Lunar Igneous rocks and the nature of the lunar interior. 

J. F, Hays and D. Walker. Proc. Sovlet-Amerlcan Conf. on 
Cosmochemlstry of the Moon and Planets, Moscow. The Lunar 
Science Institute (in press), 

2. Origin of tltaniferous lunar basalts. D. Walker, J. Longhi, 

T. L. Grove, E. M. Stolper, and J. P. Hays. Geochim. Cosmo- 
chim. Acta (in press) . 

3. Differentiation of a very thick magma body with Implications 
for the source regions of mare basalts. D. Walker, J. Longhi, 
and J. P. Hays, Proc. Sixth Lunar Scl. Conf., Geochim. 
Cosmochlm. Acta, Suppl. 6, Pergamon Press (in press). 

4. Crystallization history of lunar picrltic basalt 12002: phase 

equilibrium and cooling-rate studies. D. Walker, R. J. 
Kirkpatrick, J. Longhi, and J. P. Hays. Geol. Soc. Amer. 

Bull, (in press) . 

5. Direct determination of the quartz-coeslte transition by In- 

sltu x-ray measurements: a discussion. J.P. Hays. Contrlb. 

Min. Petr, (in press) . 
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Origin of titaniferous lunar basalts 


David Walkuk, John Longiii, HdvVard M, Stolpi r, 

Timothv L. Grovi; and Jam1:s F. Hays 

Center for Eiirih and Planetary Physics, Hoffman Laboratory, Harvard University, Cambridge. MA 021 l'.S.A. 


(Receii'ed 13 Aiiyust 1974; accopted in rcvisiul fonn 31 Oixeinhi'i- 1974) 

Abstract Delineation of low pressure phase equilibria in the composition space relevtim to titaniferous 
lunar basalts demonstrates a significant degree of control by those cquilibriti on the compositions 
of the ba.salts. The existence of two distinct chemical groups of basalts (high and low-K) which <ainnot 
be related one to the other by fractional crystalli/.ation at any prc.ssure, suggests that melting is respon- 
sible for the two groups. Consideration of the pressure shift required to produce the differencc.s between 
groups constrains magma segregation to have occurred in the outer 150 km of the Moon. It is dillicull 
to relate low-Ti and high-Ti basalts to the same source region. The preferred source region of high-Ti 
ba.salts, based on phase equilibrium considcrtitions. Is a late ilmenite-rieh cumulate produced from 
the residual liquid of the pri-mordial differentiation of the outer portions of the Moon. This ilmcnite-nch 
layer is .sandwiched between the lunar feldspathic crust and a complementary mafic cumulate. 


ISTROOI.TTION 

Thb KlRST Apollo mission to the Sea of Tranquillity 
returned a suite of highly titaniferous basaltic rocks 
and mechanical degradation products of such rocks. 
Within six months a number of contheting hypotheses 
for the origin and significance of these rocks had 
appeared in print. After four other Apollo missions 
had revealed a rich variety of lunar rock types and 
established a framework of geochemical and geo- 
physical constraints on lunar structure and history, 
the Apollo 17 mission again returned titanium-rich 
lunar basalts resembling those of Apollo 11, but more 
mafic in character and texturally more diverse. In 
view of the wealth of new information now available 
and our improved understanding of lunar structure 
and history, it seems useful to re-cxaminc the ques- 
tions posed by lunar titaniferous basalts and to 
attempt to resolve some of the points of controversy 
that have surrounded these intriguing rocks. 

These rocks have been universally recognized as 
volcanic igneous rocks that flooded low-lying por- 
tions of the lunar .surface about; 3*7 gy ago. Leading 
hypotheses for the origin of the magmas that pro- 
duced the rocks are: 

( 1 ) direct partial melt.s of the primitive lunar inter- 
ior at depths of 200 400 km. and, as such, capable 
of revealing the mineralogy and chemistry of their 
source region (R I NGWOOD and Essi-;!sii, 1970). 

(2) residual liquids of advanced stages of fractional 
crystallization near the lunar surface, probably in 
gigantic lava lakes (O’Hara et aU 1970). 

(3) partial melts of evolved regions of the lunar in- 
terior modified to greater or lesser degrees by near- 


surface fractional crystallization and/or volatile losses. 
Various depths of origin and a variety of constraints 
on the source materials have been proposed (Pmi- 
1 ‘OTTS and ScHNiiTZLiiR. 1970; Smith et al., 1970; 
W(X)D et al., 1970; Bri:tt, 1973; Prinz et al, 1973; 
Taylor and JaKi^ 1974tv. 1974 /j; Wakita el uL 1974; 
Walkhr er fl/., 1973). 

We will attempt to .show here that hypotheses (1) 
and ( 2 ) are extremely unlikely to be correct and that 
one or more variants of hypothesis (3) arc consistent 
with known facts. 

t’iii .Misnn AND PiiASi i;(;!i ii iimiA 

The titaniferous bastilts form a distinctive lunar 
type easily distinguished from other lunar magmatic 
types on the basi.s of TiOi in cxcc.s.s of fi wi and 
abundant (10 20 per cent) modal ilmcnitc (Paimki 
and Bi NCT-, 1973). Several subdivisions are recogniz- 
able on the basis of chemistry and correlative minera- 
logy. JAMI-.S and Jackso.s (1970) have .summarized 
relations from Apollo 11 where two major subgroups 
were noted: an ophilic-textured group which was 
chemically richer in ALO 3 and poorer in TiOj, 
and REE than an intersertal-texturcd group. The.sc 
are the low-K and high-K groups, respectively, iii the 
present discussion. The Apollo 17 mission sampled 
more Ti basalts, chemically similar to the Apollo f l 
low-K group but with some more malic members 
which also had TiO^ contents approaching the high- 
K group, although their chemistry did not overlap 
in other respects (LSPET, 1973). These samples exhi- 
bited a wide range of textures in contrast to the 
Apollo 11 low-K ophltic suite. The Apollo 17 Ti 
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busuits arc considered part of the low-K suite of 
titaniferous lunar basalts for the present discussion. 

Experimental studies of phase-equilibrium at low 
pressures in Apollo II compositions have been 
reported by several laboratories, These studies estab- 
lished some facts that bear on the origin of the rocks 
concerned: 

(1) the observed mineralogy and texture of the 
rocks is consistent with crystallization at low pres- 
sures and very low fugacities of oxygen and H 2 O 
(Muan and Schairhr, 1971 ; Smith et ai, 1970; Tuth- 
ILL and Sato. 1970; Weill et ai , 1970), 

(2) Rocks of the low-potassium or ophitic group 
were close to multiple saturation with olivine, pyrox- 
ene, playiociase , and/or armalcolite at the liquidus 
(Biggar et al ., 1971; O’Hara et al .. 1970, \974 a . 
imb ). 

(3) Rocks of the high-potassium or intersertal 
group had olivine, pyroxene or armalcolite as a 
liquidus phase, but plagioclase appeared relatively 
late in the crystallization sequence (Ringw(X)o and 
Essene, 1970; Ring W(K)d and Green, 1972). 

(4) Soil compositions matched neither group of 
basaltic rocks nor any mixture of them. Soil composi- 
tions had plagioclase with or without olivine or pyr- 
oxene on the liquidus (O’Hara et aL 1970, i974u), 

(5) No simple scheme of low-pressure fractional 
crystallization could produce high-potassium liquids 
from a low-potassium parent or vice versa. Complex 
schemes are possible (O’Hara et n/„ l974o). 

In addition to these experimental constraints, a 
number of geochemical constraints were established; 

(6) high-potassium rocks are enriched in elements 
such as K, Rb, Cs, REE by factors of 2 5 relative 
to low-potassium rocks (e.g, Tera et n/„ 1970; Cast 
et ai , 1970). 

(7) Rare earth elements show near-chondrjtic rela- 
tive abundance patterns except for pronounced nega- 
tive europium anomalies with absolute abundances 
levels 40-100 times chondritic (e.g. Philpotts and 
SCHNETZLER, 1970; Haskin et ai, 1970). 

(8) Rb/Sr model ages are 4-6 gy for low-potassium 
rocks and-3-7 gy for high-potassium rocks (Papanas- 
TAssiou and Wasserburg, 1971). ' 

While it was clear that the mineralogy of these sam- 
ples was the result of near-surface crysUillization, it 
was not clear to what extent the chemistry of these 
rocks was a res??;t of near-surface fractional crystalli- 
zatiort. The inability of the above observations to con- 
strain the problem sutlicicntly was partly a result of 
the narrow range of major clement chemistry in the 
Apollo 11 suite. The return of the Apollo 17 samples 
has somewhat extended the compositional range of 
the high-TI suite. 


In order to document possible near surface frac- 
tional crystallization, which might affect the composi- 
tions of the titaniferous basalts, we have performed 
melting experiments in iron capsules in evacuated 
silica tubes on crystalline powders of Apollo 17 
basalts 7(K)17 and 70215 (Longhi et ai, 1974), Micro- 
probe analy.sis of quenched glasses, coexisting with 
cry.stals, allows delineation of saturation curves in the 



oxeiic' is used to discuss the tituniferous lunar basalts. For- 
rnulac for locating molar oxide analyses in this Ictralicdron 
have been given in LoNOiii ci ui (1974, Fig. 3). Many coni- 
posiiion.s of interest have negative FmO vaUie.s, Projection 
of one sucli point is illustrated with an arhiirary symbol 
(an arrow). FmO =; (FeO -t- MgO + CaO + AUOj SiOj 
TiOzl/lFeO + MgO -t- AljOj) In molar unius. This dificr- 
entialioii index has been plotted against Fe/Fc + Mg in 
the lower portion of the figure to allow separate inspection 

01 the variation of this parameter. Experimental liquids 
from LoNcnii cr «/. (1974, Table 2) appear as dots along 
the dashed curve of this lower inset. The temperature range 
of these cxpcrimcnlul liquids, which arcahso shown in Figs. 

2 4, is 1197 1125 C. Decreasing FmO corresponds to de- 
creasing temperature in all ligUi'Cs, 

Compositions of Apollo 11 and 17 materials in this 
(Igure and Figs. 2 5 arc taken from AciRH.i t'/ (i/. (1970), 
Comi’STon h/. (1970), Ma\wiu. cf <d. (1970), RiNdWooiv 
and Essi'N (1970), Warn'i R (1971). Di nvas c't al. (1974). 
Nava (1974). Rtiooi s e/ a/. (1974), and LoNfiHl vi ai (1974). 




Origin of litanifcrons lunar basalts 


3 


portion of composition space relevant to these mater- 
ials at low pressure. A graphical representation of the 
results requires projection since there is more com- 
positional variation present than can be completely 
described by three or four components, A projection 
scheme which can be used to compare the composi- 
tions of the basalts with the curves of multiple .satu- 
ration has been presented by Walki:r ef «/. (1974) 
and Longhi et at. (1974). It is illustrated in Fig. 1. 
Details of the expcrimi^ntal methods, the composi- 
tions of experimental pblses, and an analysis of the 
strengths and shortcomings of the projection system 
can be found in Lonc-h> et at. (1974). Molar units 
are used throughout. FeO and MgO are combined to 
form FmO, but Fe/Fe -1- Mg variation is specifically 
considered in the inset included in Fig. !. FmO, TiO,, 
CaO, AI3O3, and SiOn are recalculated Into ‘minerals’ 
FmO, FmTiOj {‘ilmenite’), CaAUSiaOg {‘an’), and 
CaFmSijOfi plus FmSiOa (combined as ‘pyroxene’). 
These ‘minerals’ are then used as coordinates de- 
scribed in the captions to the various subprojections. 


Neglecting K.O and Na^O in the projection system 
has the desirable effect of suppressing compositional 
variations introduced by alkali volatilization; how- 
ever, care must be taken to note that all feldspars 
do not plot at CaAl2Si20(, as a consequence. 

Crystallization behavior of 7(K)I7 and 70215 as a 
function of pressure has been presented by Longhi et 
at. ( ! 974), The resu!t.s of this work which are relevant to 
the present discussion are that these compositions arc 
in liijiiid^s equilibrium with olivine, low-Ca cHnopyr- 
oxcneij^m'cnile,ancbqsmd in the pressure range 5-7*5 
kbar, Less olivine-normative materials such as the 
Apollo 1 1 high-K suite arc in liquidus equilibrium 
with this assemblage at pressures less than 5 kbar 
(O’Hara et at., 1970; WAbKijR, 1972). These pressures 
up to 7-5 kbar correspond to depths of up to 150 
km in the Moon. Further experimental results 'con- 
firming the applicability of these diagrams to lunar 
basalt compositions arc presented in the Appendix. 

Detailed inspection of the low pressure phase rela- 
tions for evidence of fractional crystallization control 



Fig. 2. The image plane of thi.s .subprojection is the front face of the tetrahedron of Fig. I. The 
points between the heavy dashes arc the experimental liquid compositions reported in Lonoiu vi 
«/, (1974, Table 2). The short-da.sh line between armaleoliie and ilmenite indicates the po.sition of 
this reaction on the trace of the oxide saturation surface and does not represent a tine scparuiing 
those pha.se volumes in (he oxide-over.suiurated region of the composition space. The dashed lines 
separating the olivine and pyroxene volumes are drawn on the basis of experiments on 1 2(K)2, 15065, 
and 15555 (CJROvt- el ai, 1973; Longhi ei uI.. in preparation). Dashed and dotted lines emanating 
from olivine trace oiivine fractionation paths, Low-Ti bulk compositions arc taken from C'liAt'i'LLL 
and Gri un (1973), Comi’Ston ei a/, (1971), CuTTirtA cv a/. (1971, 1973), Kosiimo and Hahami ra 
(1971), Rhoiji:S and Huu»ari> (1973), Sriru; and S.M1TH (1972). Stoi.ri-r {|9'A?), and Wii.Lis et a/. 
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on the compositional variation of the ivianirerous 
basalts is necessary before a discussion of the high 
pressure phase relations can assume any genetic sig- 
nificance, Such as inspection foltows. 

DtSCl'SSION OH PllASl; EOlilUDKIA 
ANO CUIMlCAl, VARIATION! 

One possibility to evaluate is that variability is 
caused by addition or subtraction of a single crystal- 
line species within its liquidus volume. Longhi t>( ai, 
(1974) considered this problem and concluded that 
only minor aspects of the compositional variation 
shown could be described in terms of any single crys- 
talline species. Figures land 3 show olivine fractiona- 
tion paths through the orange soil composition 74220, 
and indicate cither that equilibrium crystallization of 
olivine from 74220 or accumulation of ~Fo <70 in 
olivine-siituratcd liquids similar to some but not all 
parts of the Ti basalt suite might have occurred (see 
Appendix). Roedder and Weiblhn (1973) have pre- 
sented evidence that the 74220 composition existed 



Fig. 3 . The image plane of this projection is the rear face 
of the tetrahedron of Fig. 1. Dots along the saturation 
boundaries arc experimental liquids from Longiu ct ai 
(1974, Table 2). The composition plane can be divided into 
regions according to the silicate {plagioclasc. olivine and/or 
pyroxene) which coexi.sts with armalcolite and/or ilmcniie. 
The dotted line separating armalcolite from ilmenite satu- 
ration regions is largely schematic except in the vicinity 
of the actual experimental liquids. This projection behaves 
very much like a pseudoternary liquidus diagram for these 
materials for the physical conditions of these experiments. 
It is useful in delineating crystallization scqucncc.s' and c.sti- 
mating the proximity of saturation boundaries. It .should 
only be applied with caution to other materials, 


as a liquid with suspended microplienocrysts of 
~Fo^j( 0 . This suggests that the former possibility 
would be the more likely. However. LSPET (1973) 
on the basis of Rb, Sr, Zn, and Cl abundances have 
.sliown that neither process can in fact Ik responsible, 
since these elements are noticeably enriched in 74220 
relative to the ha,salts rather than depleted, as. would 
be expected. 

On the other hand, Lo^^-al| ct ai (19H Fig. 2) 
have shown that there is continucits (XHmpositional 
trend from Apollo 17 to Apollo II tow-K basjill 
which can be described to a very good approximation 
in terms of the multiply-saturated liquid line of de- 
.scent ob-served for Apollo i 7 experiments. Since this 
linc-of-descent was produced in a mild vacuum, there 
appears to be evidence for crystal liquid differentia- 
tion processes acting at a similar low pressure, pre- 
sumably near the lunar .surface. It should be noted 
that this trend does not include the Apollo 1 1 high-K 
suite. Projecting from alkali feldspar components in- 
stead of NaiO and K 2 O does not substantially affect 
the relative dispositions of these composition 
groups. 

A subprojection is seen in Fig. 2. The experimental 
liquids are saturated with armalcolite and/or ilmenite 
and .so ilii.s heavy-dash curve represents the Ti-oxide 
.saturation boundary in this portion of composition 
space. The volumes for olivine and/or pyroxene crys- 
tallization are indiaited on the Ti-oxide undersatur- 
ated side of this curve. The relative positions of these 
volumes have been determined by work on Apollo 
12 and 15 iow-TI basalts (Grove et «/., 1973; Longhi 
et ai, in preparation). Itcan be seen tea lirst approxi- 
mation that the Apollo 1 1 and 17 basalt compo.sition.s 
cluster along the curve of Ti-oxide .saturation. The 
most noticeable deviations arc seen for the Apollo 
11 high-K suite. The deviations arc such that these 
materials are over-saturated with respect to Ti-oxide. 
This compositional distribution contrasts strongly 
with that of the low-Ti mare basalts also dusplaycd 
in Fig. 2, which are undersaturated with respect to 
Tl oxides. Olivine fractionation paths arc indicated 
with dotted and/or dashed lines. Extreme fractional 
crystallization of low-Ti liquids involving olivine then 
pyroxene (and plagiodase) is capable of producing 
Ti-oxide-saturated liquids, A liquid produced in the 
experimental crystalli/aiion of 12002 (Grove vt w/., 
i973) is indicated by the diamond at the end of the 
arrow. These residual liquids bccorne timenite-satu- 
ralcd but only at temperatures below li25C and 
below the armalcolite .stability held even for the most 
tiUiniferous of the iow-Tl basiilts, The faet that armal- 
colitc- (not ilmenite )-.siitura led compositions such as 
70215 appear to have existed in the liquid state at 
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~ 1200'^C very sirongly suggests that such liquids are 
not the dregs of fractional crystallization of these or 
similar low-Ti magmas. Furthermore, fractional crys- 
tallization of low-Ti basaltic liquids will produce Ti- 
oxide-saturated liquids with high Fe/Fc ■+ Mg (>() 8) 
unless there is simultaneous separation of a metal 
phase, If meia! does not .separate, the liquid pro- 
duced is too iron-rich. If metal does separate, olivine 
becomes unstable as silica activity increases. The 
suggestion of O’Hara X'/ a/. (1970) that titaniferous 
basalts represent the residual liquids of advanced frac- 
tional crystallization cannot be sustained unless the 
parent is already highly titaniferous so as to produce 
armalcolite siituration. Low-Ti basalt is not a satisfac- 
tory parent. 


C)n the other hand, the hypothetical parent liquid 
must not be so titaniferous that major precipitation 
of Fc Ti oxides occurs with consequent depletion in 
Cr (Ri.\gw(K)d, 1970). The range of plausible parent 
liquid compositions is therefore narrov,^ly restricted. 
Among the returned lunar samples only the “orange 
soil' (74220) falls within this range and it is unsatisfac- 
tory on other grounds (Appendix; LSPET. 1973). We 
conclude that the titaniferous basalts are unlikely to 
be direct products of extreme fractional crystalliza- 
tion. 

The approximation that the titaniferous bastilt 
compositions arc saturated or oversaturated with 
titaniferous oxide makes possible another useful pro- 
jection, Figure 3 may be separated into silicate satu- 



Fig. 4. Projection from FmTiO.t on to plane SiO^ An Px, This contains the sjime information as 
tlic projection of Fig, 3, but it is recast into a more familiar form. The Px and An molecules arc 
the same but SiOa is used instead of FmO: the ternary format uses “olivine’ (iiot ideal because of 
CaFmSijOf, component in the pyroxene) as a coordinate. Molar oxide analyses may be plotted in 
this projection according to the following formulae; 



(SiOj Gaa AljOj + iTiO ] FcO iMgO) 
“ SiOr-GiO 


CaAl2Si20g 


M2O3 
Siol CaO 


‘Olivine’ 


FeO -i- Mga TiOi 
2 X bSiO, CaO) 


Note that FmO is now at infinity below ‘olivine’ on the 'olivine* silica join. For comparison the 
saturation boundaries for other systems are shown when recalcuialed in this projection. There is u 
systematic shift of piagioclase-fcrrqmagnesian silicate saturation curves away from plagloclasc with 

increasing Fc/Fc + Mg of the systein. 
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ration fk*ld& which cxicxiRt with an oxiiic (ilnwnib; 
at>d or arnukolitc). The projection is from FmTiOj 
on to the KmO pyroxene CaAIjSijOi, plane and 
hehavex like a ternary liquidus diagram in that crys- 
talli/ation paths and the order of phase appearance 
may be deduced. New experimental evidence which 
confirms the applicability of these boundaries to high- 
K and other materials may be found in the Appendix. 
(The uimpositions of the basalts and the saturation 
boundaries are recast in a slightly more familiar ter- 
nary format in Fig 4.) 

Several impi>rtant conclusions can be drawn from 
consideration of the positions of the titaniferous 
basalt compiisitions relative to the saturatKin boun- 
daries in Figs 2 4. The compositions seem to be dis- 
posed along the surfaces of multiple saturation. This 
observation supports the contention of O'Hara et al. 
(19701 that the Apollo II basalts were ‘near-cotectic’ 
in compiisition. and indicates that these compositions 
are controlled to a significant degree by crystal liquai 
fractioiation at low pres.surc. Note that ‘multiply- 
*saturated liquids' arc not synonymous with ‘plagio- 
clase-saturated liquids' although plagioclase-saturated 
representatives of the titaniferous basalts are found 
(see Appendix. 750.^5). Note also that our detailed 
confirmation of the experimental work of O'Hara et 
ii/. (1970, 19741 and Biguar et al. (1971. 1972) is not 
to be construed as an endorsement of their conclusion 
that advanced fractional crystallization near the lunar 
surface is the crystal liquid equilibrium proc'css re- 
s|x>nsible for the compositions of titaniferous basalts. 
It has been argued above that an origin by advanced 
fractional crystallization encounters difficulty in 
explaining the presence of magmata; armalcolite. high 
Cr and V (Kingwikii). 1970) and rebtively modest 
Fe Fe Mg. unless rather severe restrictions arc 
placed on the parental liquid. We will argue that 
melting priKcsses are a more reasonable means of 
gerterating the bulk chemistry of two distinct groups 
of high-Ti basalts. This is to say that the 
demonstration of low pressure crystal liquid frac- 
tionation control is not sufficient to den'K>nstratc that 
all the control of chemical variation is by fractional 
crystallization. We (k\ however, recogni/c that these 
compositions are controlled by low pressure equili- 
bria to a certain extent and that fractional crystalliza- 
tion may interrelate various types within the titani- 
ferous suite, since only relatively minor fractionation 
IS required. For instance, it is possible to consider 
70215 a rapidly chilled sample of a st^rt of magma 
which may have cooled and differentiated by olivine 
and armalcolite fractionation to give 750.f5 (the 
Apollo 17 basalt ptxjrest in FmO). Brown et ai 
(1974) have suggested such a petrologic relation may 


exist between various numbers of the Api>llo 17 suite. 
A similar case can be made relating the members of 
the ApeMo II low-K suite except that quickly chilled 
samples of likely parent iH^uids were not sampled 
Figure 5 sht>ws minor and trace element data compat- 
ible with this interpretation. The umiptfsittonal wna~ 
tkni of Apolht 17 anJ Apitllo II /<m-A titaniferous 
basalts is eontinmms, otvrlappom. and expliciihle as the 
result of mtniest anumnts of fractional crystallizatuMi 
of (divine and armalcolite folhm ed by ilnhmite then pla- 


HTAl 



Fig. 5. Minor and trace element data for Apollo 1 1 basalts 
are plotted against the FmO value used as ordin4ilc in 
the inKt to Fig. I. This parameter serves as a differentia- 
tion index with increasing differentiation to the right The 
ophitic-intersertal classification follows Jami s and Jai xsos 
(1970). Within the ophitic suite there are crude trends to- 
ward decreasing Cr^Oj and increasing KjO. La. and Sm 
with differentiation These trends are consistent with frac- 
tional crystallization of armalcolite. ilmenite. olivine, 
and plagKKiasc. No such trends are evident in tfw interser- 
tal suite, suggesting that SiO^ analytical problems may be 
artificially inflating the differentiation index of tlw three 
samples to the right (sec test). Symbols are average values 
reported in Warnik (1971). Bars indicate the extreme 
values reported. Oita density is not uniform 
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tfHHluH- (.V/ thal llh' Jerh'uliiv arc plmiMh 

du.st'-saiuruh'J). Howrivr,Jttr //k* rcuswis milluH'J ahtnr 
it is unlikt-ly that the parental lUpad is itself tlx- result 
of extensile fraetionai crystallization. 

It appears unlikely that the Apollo 1 1 high-K suite 
can be simply related by low pressure fractional crys- 
tallization to either the Apollo 17 or Apollo 1 1 low-K 
suite. This point has been made successfully by Jamks 
and Wright (1972) by consideration of the crystalli- 
zation sequence within each group and it can be seen 
in the projections of Fig. 3 and 4 that the high-K 
suite is not plagioclase-saturated as it should be if 
It were the fractional cTystalli/ation product of pla- 
gioclase saturated Apollo 1 1 low-K liquids. There is 
a curious aspect to the compositional variation of the 
Apollo 1 1 high-K suite; it appears to parallel the 
variation of pyroxene plagicK'lasc ilmenite-saturated 
liquids If this distribution is in fact caused by satu- 
ration with these phases, one most postulate some 
set of conditions w here the equilibria (stable or meta- 
stable) arc shifted in favor of less feldspathic liquids. 
Pressure increases or hydration would shift the rela- 
tions the wrong way. Alkali devolatilization would 
shift the relations the correct way but it w ill be argued 
later that the magnitude of the effect is not sufficient 
and. furthermore, the high-K suite would appear to 
have sustained less alkali loss than tlie low-K suite. 
Alternatively this distribution may be an artifact of 
SiOj analytical unc'ertainty in the analyses gleaned 
from several sources. The three low FmO (high SiOj) 
points in the distribution are data obtained by acti- 
vation analysis and spark source mass .spectrometry, 
in contrast to XRF and wet chemical determinations 
on which the rest of the distribution is based. With 
reference to Fig. 5. the uniformity of minor and trace 
element abundances within this group (as opposed 
to the rudimentary trends within the low-K suite 
suggestive of crystal fractionation) makes it more 
likely that this trend is the result of small SiOj ana- 
lytical uncertainties. If this is the case it is not necess- 
ary to postulate shifted plagicK'lase pyroxene equili- 
bria to explain an artificial distributK>n. but it is then 
necessary to postulate a separate parental liquid for 
the high-K suite. This separate parental liquid cannot 
Itself be the product of advaiKcd low pressure frac- 
tional crystallization for the sanne reasons discus.sed 
for the parent to the low-K suite. 

AIXAII ViH Atll I/Anos. PlAiiMNlAV Sahhatmis. 

AM) ‘AVIRAta' MaOMA ('oMIHIfUTHIN 

Much has been written concerning the possible role 
of plagUK'lase fractionation in the production of the 
titanilcrous basalts, and much of the controversy has 


been over the location of the plagRtclasc saturation 
surface. It can be seen in tlie projations of Figs. 2 
and 4 that the plagHiclase saturation surlace 
encroaches upon sonw of these compositions. This 
circumstance has been repeatedly emphasized by 
O’Hara et al. (1970. |974<i) and Hkr.ar et al. (1972) 
and IS not open to further questR>n. It can also be 
seen that piirts of the low and high-K distribution 
are not immediately plagioclase-saturated. These 
compositions are the ones we consider to be parental 
to the rest of the low and high-K suites through 
limited fractional crystallization. 

RimiW(ioi) and CiRbrs ( 1972) (and others, including 
O'Hara et aL 1970) sluiwed that high-K composi- 
tions (e.g I(N)I7. 10072) did not have plagioclase on 
the liquidus. An extended controversy developed with 
one group emphasizing the plagKK'lasc-saturated 
compositions and the other group emphasizing the 
plagioclasc-undersaturated compositions. 

O'Hara et al. (1970) made the suggestion that the 
undersaturated materials could be made plaguK'lase- 
saturated by restoring alkalis and other volatile con- 
stituents presumably lost by volatilization at the lunar 
surface. Ringwood and Grttn (1972) performed ex- 
periments to test this contention and found that addi- 
tion of alkalies suppresses .atpprrggy plagitKiase 
appearance and makes plagioclase saturation more 
remote. We have performed alkali addition exper- 
iments on 70215 and obtained the same result. The 
pyroxene and plagioclase saturation curve nK>ves 
closer to plagioclase and hence plaguK’lase saturation 
is more remote from the compositions of interest (sex* 
Appendix). If addition of alkali makes plagiiKlase 
saturation more remote, then loss of alkali should 
increase the proximity to plagioclase saturation (see 
Fig. 6). 

It is worth pointing out that the two microprobe 
analyses of Biggar et al. (1971) for Apollo 1 1 exper- 
imental residual liquids are in exact agreement with 
the saturation curves presented here. Bkigar et al. 
(1971) did not perform their experiments in vacuum, 
as we did. so their liquids have retained roughly three 
times the alkali content of those reported here. This 
implies that the magnitude of the alkali-loss effect is 
not sufficiently large to aecount for the deviations 
from cotectic character since experiments in whK:h 
substantial alkali loss has been sustained, relative to 
a low starting value, do not move the curves percepti- 
bly. The effect denwnstrated by RimiW'irid and 
Grii n (1972) was only accomplished by substantial 
addition (2-5 per cent) of alkali. 

A related issue has been raised by O'Hara et al. 
(I974u). Their contention is that tlK hand specimen 
compositions are not truly representative material but 
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Fig. 6. FxpcrimcnUil rcsult> arc displayed on the proicclion of Fig 4 Bulk compositions of the starting 
materials and the experimental crystallization sequences may he found in Tables I and 2. The shift 
of the siituration houndarK's as a result of changing physical conditKins is indicated Addition of 
alkalies shifts the pyroxene plaguKlasc hound iry towards plagiiK'lase. Compositions of alkali-enhaiKcd 
liquids indicated by diamond symbols are four^J in Table I. Increasing pressure shifts the olivine pyrox- 
ene K>undary towards nusre oiivine-normativc compositions. 


rather that statistical studies of lithic fragment com- 
positions in the sr>ils and breccias are more mean- 
ingful prov ided that the effects of contamination with 
highland debris can be removed. The ‘average mag- 
nvi' they obtain is indeed plagioclase-saturated. This 
approach renders arguments over the hand specimen 
chemistry, and whether it has been devolatilized, irre- 
levant. since the shift from plagioclase saturation is 
not brought about by a change in the saturation 
'^>i*V)undary but by a change in hand specimen composi- 
tion. i.c. by partial accumulation of small, pet^ra- 
phically cryptic, amounts of ferromagnesian silicate 
and oxides. This is to say the hand specimen composi- 
tions do not represent liquid compositions. 

In contrast to this approach, we and most others 
have made the tacit or explicit assumption that the 
hand specimens do represent various lUfuui composi- 
tions since there is no obvious petrographic evidence 
to the contrary. Indeed several Apiillo 17 samples (e g. 
70215 and 71.^69) show quench textures, indicating 
cooling from a totally liquid state. The inference that 
the liquid had the siime composition us the hand 
specimen is insecure, however, as fuis been pointed 
out by O'Hara (personal communication). Mechani- 
cal transfer of cryptophemKrysts or ditTusioiuil 


transfer through the liquid might substantially alter 
the solid a>mposition if the transfer 'ate exceeds the 
cooling rate. If this be the case, it becomes necessary 
to invoke some special nwchanism to explain the 
compositional trends show n by both high-K and low- 
K hand-specimen groups and their relationship to the 
experimentally determined cotectic curves. Wc prefer 
to maintain the assumptam that compositions of 
hand-specimens showing glassy or queiK'h textures 
arc the most reliable guide to magm;i compositions. 

If we have correctly identified the low-K parental 
compositions, limited fractional crystallization of 
these can potentially produev a large volume of resi- 
dual cotectic magma. That some has been produc'ed 
is clear from samples such as 75().^5 (see Appendix). 
This derivative, possibly plagioclase-phyric. material 
could form up to X() per cent of the crxiling unit if 
differentiation were ctficient. (Tvidcntly the differen- 
tiation was not completely ett'icient since most of the 
range of intermediate compositions between parent 
and cotectic have been found in the hand specimen 
suite.) O'Hara's paradoxical discrepancy between 
hand specimen compositions and compositions of 
lithic fragments in siiils and breccias suggests the 
existence of some sampling bias. Since the cooling 
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unit thickncsik and sampling mechanics arc unknown. 
It appears equally probable to us that the sr>il could 
be skewed toward the derivative compositions instead 
of toward the parental compositions if the cooling 
unit contains subequal parts of preserved parent and 
derivative cotectic material. 


OkKJIN OI- PARI NTAt LKyillDS 

We have shown here that the titaniferous lunar 
basalts exhibit compositional variations within the 
suite that appear to be controlled by low-pressure 
crystal liquid equilibria. We have also argued that, 
although the relationships among the samples in the 
Apollo 1 1/17 low-K suite appear to be controlled by 
modest anu^unts of fractional crystallization at or 
near the lunar surface, the parental liquids themselves 
arc not likely to be products of advanced low-pres- 
sure fractional crystallization. Kurthermi^re. there 
appears to be no plausible scheme of fractional crys- 
tallization capable of relating the high-K Apollo 1 1 
to the low-K Aptillo 11/17 suite. 

We arc therefore forced to seek an origin for at 
least two parental magmiis by partial melting some- 
where within the lunar interior. 

RiNtiwooi) and f sM si (1970) attempted to locate 
the depth of origin and source material for the Apollo 
1 1 suite by means of high pressure studies of a ‘model 
basalt comptisition’ based on an average of Apollo 
II high-K and low-K rock analyses. Unfortunately, 
this miHlel composition was prepared at an early date 
when only preliminary emission spec trograp hie data 
on rock compositions were available, and as a conse- 
quence its MgO KeO ratio is higher that that of aver- 
age Ap<dlo 1 1 magma, although w ithin the range of 
values reported (Kingwood. personal communication). 
In most other respects it lies near the pyroxene pla- 
gKKlase ilmenite-saturated experimental liquids (Figs 
I .^) and far from the comfx^sitions that we would 
regard as possibly parental and least likely to have 
been affected by low-pressure fractional crystalliza- 
tion. We therefore reject the main conclusion of 
Kinow(x)I) and Fs.si m (1970) and RiMiWixiD (1970). 
which was based on the study of their model basalt, 
that titaniferous magmiis arc produced by partial 
melting of primitive lunar material at a depth of 200 
400 km. even though nu^st of the general arguments 
presented in these papers are valid 

We have reported elsewhere (LoMiiii et uL 1974) 
on high-pressure studies of 702 1.S. a composition 
which we regard as a possible parent magma for the 
low-K suite. Such a liquid would be in equilibrium 
with a crystalline residue of olivine, pyroxene. 


^ilmenite .iiiJ gMimibl!> 1 1 i > <prw»l at some depth near 
120 km. A depth of origin greater than l.^f km would 
suggest unrealistically large degrees of partial melting 
A less olivine-rich high-K parental liquid would be 
m equilibrium with a similar residual assemblage at 
even shallower depths. Such an episode of magma 
segregation at modest depth precludes the possibilit) 
of deducing the nature of the deep lunar interior on 
the basis of high pressure study of the titaniferous 
basalts. 

An ilnK'nite-rich assemblage is unlikely to represent 
primitive lunar material, and the shallow interred 
depth of origin strongly suggests tivit the source 
material has been subjected to the early magmatic 
differentiation that produced the feldspathic lunar 
crust. A piissible primitive liquid .such as that sug- 
gested by W'alkik el ul. (I97.f; Fig. .V pt. A) with 
0^9 wt TiOt and Fe (Fe Mg) (f24 (molar) 
would begin to crystallize ilmenite when approxi- 
mately 90 per cent crystallized (reversible path). The 
residual liquid at this stage would fuive Fe (Fe 4- 
Mg) = 0-45 (F50 and would be in equilibrium with 
olivine (Fo,,,, 77 k pyroxene, ilmenite and plagioclase 
(Departure from reversible crystallization would yield 
more iron-rich residual liquids.) Such a liquid would 
produev sinking cumulates of olivine, pyroxene, and 
ilmenite. floating plagioclase. and LiL enriched resi- 
duum. Later remelting and mixing of residuum with 
melted mafic cumulate at variable depths of up to 
150 km would yield titaniferous parental basalts cap- 
able of further near-surface fractionation to produce 
the observed Apollo 1 1 and Apollo 17 riK-ks. Variable 
proportions of residuum to cumulate phases in the 
molten products might explain some of the observed 
minor element correlations. 

Similar mixk'Is have been proposed by others (e.g 
Philwtts and StHM r/.t.i r. Smith ei uL Wood ei 
uL all 1970; Britt. 1973; Taylor and Jams. I974u. 
1974/); W AKITA I’l aL 1974) to explain the Fu-ano- 
maly and other details of minor element geiKhe- 
mistry. A vexing question has been whether or not 
plagioclase was present in the source nKitcrial. If pla- 
gioclase is not present in the post-melting residue 
then the source nuiterials must possess a strong in- 
trinsic hu-anomaly. Such an intrinsic aiH>maly seems 
a natural consequence of the prior separation of a 
feldspathic crust. In the case of the low-K titaniferous 
rocks (both Apollo II and Apollo 1 7) the existing 
Rb Sr data yield model ages near 4 6 gy. requiring 
that no significant fractionation of Rb and Sr 
occurred after that date Hence tlK proptvsed partial 
nK'lting event at 37 gy must have quantiUitively 
extracted Rb and Sr from the source material; pla- 
gKKlase cannot have remained in the residue 
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(although it may have been present before the melt- 
ing!. This is consistent with the phase equilibrium 
data for the supposed parental liquids The high-K 
tiUiniferous riK'ks of Apollo II have model-ages near 
' ^t'*'i -the crystallization ages (PAiyANASTASSKHi and Was- 


-> Ad 1 ^^ St Kilter i. I')7I). This suggests that some strontium- 
bearing phase, presurmibly plaguK'lase or calcic pyr- 
oxene. was left in the crystalline residue. Available 
phase equilibrium daUi arc not adequate to identify 
this phase conclusively but indications are that it is 
not plagK)clase. 


RllAnONSHIH TlTAMtlK«M% MaNI 

Hasai is n> 1 .o\a-Ti Mahi Hasai is 

It has been proposed by CIkiin et al (1071. 1074) 
that low-Ti m;ire basalts represent partial melting 
products of an olivine pyroxenite sourev regarn at 
depths of 200 400 km in the MiK>n. This proposal 
was suppt>rted for the case of Aprillo 12 by Grovi 
t'f al. (107.^! but the matter is not beyond dispute, 
as is clearly shown by the arguments of BKi<iAK ci 
li/. (1071. 1074). An important problem in m;ire basalt 
petrogenesis is: can the low- and high-Ti varieties be 
related to a single source by various degrees of partial 
melting and or various depths of magma segregation? 
Griiv et al. (1074) have argued tliat it is pi>ssihle 
to relate the two types at a depth of 400 5(K) km 
on the basis of studies on Apollo 17 orange soil, a 
strongly olivine-normative titaniferous composition. 
We argue (in the Appendix) that such a parent to 
the high-Ti suite cannot prtnluce both the high- and 
loW'K suites by fractional crystallization at any depth 
and SA) we are constrained to accomplish magma 
segregation from an ilmenite-rkh source in the outer 
150 km of the Mix>n. We do not endorse the proptvsi- 
tion that bi>th low- and high-Ti basalts can be derived 
at depth from the same primitive source material. We 
arc arguing for the necessity of a near-surface, differ- 
eiuiaieJ stvurc'C region for the high-Ti basalts. Dn\- 
t AN et al. (1974) and PiiiLHon^ et al. (197.^) have 
demonstrated that there is also cvideiwe in the trace 
element abundance patterns to support the notion of 
a two-stage model of titaniferous magma production. 
RiNtiVMioi) ( I974</. 1974/)), on the other hand, has sur- 
veyed the broiider aspects of the trace element abun- 
dance patterns in nvire basalts and proposed a single- 
stage model for the production of both low- and high- 
Ti basalts from the same source material. It seems 
premature to us to construct elaborate models which 
attempt to reconcile the trace element abundance pat- 
terns of the various basaltic types when the major 
element relations are in doubt. In this connection it 
IS worth noting tliat all nvxlels proposed have ditli- 


culty reconciling the trace element abundance pat- 
terns in detail. In a word, complete consistency with 
regard to major and minor elerrcnts has not bex'n 
achieved by any model. 


I HtlU IN»\S )<• Ml 1 IIMi tN Cl Ml I AH.A 

KiMiWiHii) ( |974/>) lus raised four objections to the 
generic class of models wharh melt cumulates to pro- 
duce titaniferous magmas. The nxYst serK>us of these 
objections involves the observation that source 
regions which have lud a prior history of crystal 
liquid fractRvnation should prinluce basalts with 
strong minor elenu'nt fractionations, for example, 
fractionation within tiK* KFE. The europium anomaly 
of these basalts (its the mtxicl very well, but one might 
also expect relative depletam in the heavy RFE. This 
objection is not easily met in equilibrium crystal frac- 
tionation models and it is just this sort of objection 
which has forced Riviwtxii) ( I974u. l974/>) to aban- 
don the bulk equilibrium boundary condition if this 
constraint is renaoved the objection is less serious 
It is possible, however, to construct models in which 
cumulates depleted in the light RFF with an intrinsic 
europium anomaly might give rise to partial melts 
with apparently un fractionated RFF patterns (except 
for europium) without abandoning the bulk equilib- 
rium constraint. These mi>dels constrain the cumulate 
RFF pattern and the degree of melting required 
Much progress has been made in constructing such 
mixlels by Tayior and Jake?^ (I974<i, l974/») and 
W AKITA et al. (1974). 

Another objection involves the lack of chromium 
depletion which might be expected if the source 
region is pri>duced by extensive fractional crystalliza- 
tion. O'Hara et al. (I974«) have pointed out the 
weakness of this argument, ('hrome depletion by 
spinel must exceed chrome enrK.hment by olivine and 
plagioclase for the argument to work. Mivdeling 
chrome depletion by analogy with internal fractiona- 
tion in crystallization of high-Ti basalt hand speci- 
mens (Ristiwixu>. 1974/)) is not relevant since 
ilmcnite is also fractionating there The only signifi- 
cant Cr depletion envisaged in the fractionation pro- 
ducing the Ti basalt source region results from spinel 
crystallization which must compete with C'r enrich- 
ment caused by subtraction of olivine and plagioclase. 
The competition will not deplete C'r unless spinel is 
more than about 1-3 per cTnt of the total crystal 
extract. C'r is, of course, depleted, possibly from a 
high starting value, when the residual liquid begins 
to crystallize ilmenite But Cr is recovered in melting 
that ilmenite to produce the Ti basalts. 
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A third ohjcction regards the absence of a suitable 
heat si^urce to melt shiilkm cumulate layers alter the 
outer parts of the Moon have cooled by conduction 
If it is assumed that the radioactivity of the residual 
liquids from which the cumulates formed was sutli- 
cient to remelt the cumulate at 3-7 gy then it is hard 
to understand how the liquid could have crystallized 
in the first pluc'e. John Wcxxl (personal communica- 
tion) has suggested that a transition from a convective 
mode of heat loss through the lunar crust to a uvn- 
ductive mixk' ntiy have accompanied the dc'cline in 
impact activity at about 4 gy. and thtit this transition 
may have initiated a subsequent rise in temperature 
in the subcrustal regions. 

Smith (1974) and Wonm ( 1974) have suggested that 
VISCOUS tidal heating may be an important lunar heat 
source: and this mechanism is especially useful for 
heating the postulated source region of the titani- 
ferous basalts. This source region is sandwiclK'd 
between complenK'iitary feldspathic and mafic cumu- 
lates of contrasting density and mechiinical properties 
and may have trapped small segregatKins of uncon- 
gealed liquid residue as noted above. Such a layer 
would be an ideal locus for dissipation of tidal energy. 
It is not clear, however, that tidal energy available 
at the tinw of nvirc filling would have been adequate. 

We consider the heat source problem to be serious 
and unst)lved. It seems clear, however, that the equi- 
librium clKmical relations do not indicate a deep 
source region where the thermal and rheological pro- 
blems might be less severe 

A fourth objection concerns the feasibility of mo- 
ducing late stage cumulates within the Moon. 
possible that the late liquids will segregate and dif!^ 
entiate at a level hdow the feldspathic crust rather 
than be extruded.* This question is transitional to but 
clearly distinguishable from the one which asks if feld- 
spar dotation in a global magma ocean can be re- 
sponsible for the feldspathic crusUil rocks and the late 
residual liquids. An attempt was made to answer the 
former question by equilibrating a crystalline powder 
starting material synthesized to have the composition 
of Apidlo II soil I00S4 (Fromml ct uL 1970) at 
I2(N) C for 4] days. This titaniferous composition has 
plagKK'lase at the liquidus (WAt KiK. 1972) and can 
be used as a rough model for the residual liquids. 
The quenched product of this experiment (in vacuum 
in high purity iron container) was >95 per cent glass 
w ith plagKK'lase concentrated at the top of the charge 
and spinel and iron metal droplets at the bottom 
The plagioclase coarseiK'd upwards in the charge, and 
the clear implK'atKm is that plagioclase floats in .such 
a liquid The density of this liquid can be calculated 
from the model of BoTTistfA and Wmli (1970) to 


be 3-25 gm cm^ at 1200 C'. which is consistent with 
plagKK'lase flotation whK'h is observed experinKn- 
tally. Pstimates of average lunar highbnd composi- 
tKvn such as anorthositic norite have low pressure 
mineral assemblages which give densitK's -29 gm 
cm^. Late-sUigc liquids parental to the ilmcnite-rKh 
cumulates should therefore remain beneath such a 
feldspathK crust. 

Com 1 1 min 

Wc have shown that near-surfaoe fractKvnal cTystal- 
li/ation IS capable of explaining the compositional 
variations within the Apollo 1 1 low-K and Apollo 
17 sublltHir basalts, but that liquids parental to this 
suite are unlikely to be the result of advanced frac- 
tional crystallization. The Apollo 1 1 high-K suite can 
not be related by fractional crystallizatKvn to the 
Apollo II low-K and Aptvilo 17 subtliH)r basalts. 
Compositions that are parental to tin.* Aptillo 1 1 and 
17 suites can be generated by partial melting of a 
differentiated sourex* regKm (ol -i- px + 
the outer 150 km of the Moon. This model is capable 
of explaining the origin of and differences between 
the Aptvilo 1 1 high-K basalts on one hand and the 
Apollo II low-K and Apollo 17 subtliKH basalts on 
the other hand. The differences arise from varying 
degrees of melting of the stvurce region and variable 
depths of magnu segregation In this model the 
Sivurce region docs not contain plagKK'lase (or at least 
diKs not retain it in the residue after melting) and 
the differeiK'es result chiefly from different depths of 
magma segregation. However, if plagioclase is pres- 
ent. an alternate scheme in which the differences 
result from increasing degrees of melting is possible 
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Appyhiax 

I YPtKIMtMAI RiSI ITI AM) DlSt l ViMIS 

1 ilunift'rous husulls 

Wc report here further cxpermYcntal rcsulU relevant to 
the titanilerouik hasalt problem. Details of the experimental 
priKX'durcs nui> be found in Lonhgi er ul. (I974|. Table 
I contains the bulk compositions of starting materials used 
in this investigatmn and the compositions of selected 
experimentally produexY pluses. Kigure 6 shows these com- 
positions in the same projection as Fig. 4 Table 2 lists 
the conditions and results of the individual experiments. 
It can be seen from an inspectum of Table 2 and Fig 
b that the low-pressure saturation cur vs determined by 
microprobe analysis of residual glasses ir> 7(NH7 and 70215 
(Wai kir t’f ul.. 1^74, Lonomi ct uL I974t very closely 
predict the crystallization behavior of the three other 
titaniferous basalts reported here. This co istitutes a conlir- 
nvition of the utility of this projection. 

An important feature of these experiments is that high 
purity iron capsules are used inside an evacuated silica 
tube. The use of iron capsules introduces two problems. 
There is little direct petrographic evidence that these com- 
positions were iron-saturated at the liquidus; hence, they 
might be expected to react with the capsule and gam iron 
This eflect was sought and was undetectable in the original 
work reported by Walkir er at. (1974) and detailed in 
UisoMi t'l ul. (I‘>74) on 7()0|7 and 70215. In the present 


Table I Electron microprobe aiulyses 



rsojs.t-i* 

71bb9.21‘ 

10072. Z3* 

7M220.3/* 

ElOj 

VI. 6 

3B.7 

38.7 

19..' 

TlOj 

10.2 

13.0 

13.0 

9.VV 

AljOj 

0.7b 

B.6V 

7.56 

5.6V 

Cr,oJ 

.27 

.*5 

.30 

.6V 

Peo 

Vi.b 

lh.9 

21. V 

22.2 


6.K, 

7.87 

7.36 

IV. 1 

tno 

.30 

.i« 

.3? 

.3V 

CsO 

U.b 

10.3 

9.73 

7.10 

.(•jO 


.VI 

,*7 



9U.7C 

98.61 

96.8V 

99.2V 



Dli-l* 

215-SdJ. 

215->*,- 

215-5P--} 


op*** 

-class'* 

. last** 


alOj 

52. V 

V8.7 

V6.1 

vv.i 

TlOj 

1.10 

5.55 

6.93 

7.96 

AI 2 O 

3.53 

12.6 

12.0 

11. V 

Cr20’ 


.05 

.05 

.07 

F*0 

13 . V 

15.5 

Ib.V 

16.9 


26.2 

3.62 

V.35 

5.17 

In 


.26 

.28 

.29 

CaO 

2.27 

8.26 

9.63 

10. b 

.(■^0 


JuU 

V.51 

-liii 


98.90 

99,87 

l;o,2‘ 

99. “I 


* Hulk compositHvns of experiments abt>ve the liquidus: 
750.^5.bl low-K chemistry; 7156^,21 intermediate che- 
mistry: 10072.2.^ high-K chemistry; 7422fk.n orange 

s*^ 

*^Sce Table 2 for conditions of the experiments / 


work the ellect was undetecMble in 74220. 7SI.T^ and 
7l5bM For 10072 the efkxt pr«>duced a gain ol 7 per cent 
of the amount of iron present in the l^lk compositHut 
With the exceptHin of HNI7Z tfN.'n. the compositmns must 
have been effectively iron-saturated to mountain constant 
bulk composiimn This .irgument has alvi been made 
recently by Kivvin ( I974j with reference to low-Ti basalts 
from Apollo 15. For l(’072 we have reported tlie mildly 
iron-enhanced composr.Kvn in TaNe I and in Pig b The 
composition still falls in the range of the high-K suite of 
which It IS a member, and as such, its experimental crystal- 
lization si'quencx’ is instructive. 

The second problem introduced by iron capsules has 
becri noted by Lisust.i v er ul. (1974) Armalcolite will show 
a !arger stability field in iron capsules than in non-iron- 
iTiturated conditions. This efkct will tend to ntake armal- 
colite persist with ilmenitc to lower temperatures. 

Another feature of our experimental configuration is a 
large void-space sample ratio which allows charge extrac- 
tion with re-use of the capsule. The large evacuated void 
is presumably influential in the p;irtial loss of alkali experi- 
eiKed by our charges although the final disp«>sitH)n of the 
alkali IS unknown. By comparison with saturatmn-curve 
information in Biooak er ul. (I*>7lf the loss ol alkali reb- 
tive to the very low starting level present d«K's not substan- 
tially aflevt the positam of the saturation curves. Any effect 
introduced by the alkali loss would make the plagHKiase 
more calcic and stable at a higher temperature relative 
to a non-alkali^epleted situation. The coupled elTixt. titen. 
of our experimental perturbations is to increase armucolite 
stability to hiwer temperatures and pbgHKlasx* stability to 
higher temperatures than might be otherwise achieved As 
a result our experimental ‘gap' between tlte removal of 
armalcolite by reaction artd tlK saturation of the liquid 
with both plaguK'lase and pyroxene must be a minimum 
gap. Under other conditkvns it would be larger This urnsi- 
deratHin strengilurns the contention in the text that armal- 
colite-bearmg lavas cannot he erupted from plaguK'lase 
pyroxcne-cotectic magma bodu's. Independent indications 
of the truth of this proposition may be found in the petro- 
graphy of 75055 (Kridh HAl (ill and Wmi. I*#7H 75<H5 
(LoNciiii ef ul.. IV74k KNI44 and KMH7 (Jamis and Ja( u- 
siis. I*>70k which have cotectic bulk comptisitums. These 
samples do not contain armalcolite. and 750.^5 docs m>t 
crystallize armacolite experimentally as noted below 

W'e have given a petrographK desenptum of medium- 
grained ophitK Kisalt 75(H5 (LomiIII i*i ul . I*>74) in whK'h 
ilmenitc was rcx'ognized as the first crystallizing pbise. 
quK'kIy followed by plaguKlase and pyroxene. I xpcrimen- 
tally. pyroxene and plaguK'lase crystallize simultaneously 
immediately after ilmenite saturatu>n. The maximum tem- 
perature interval between the lujuidus and thrcv-fold satu- 
ratum (ilmenitc. pyroxene, and plaguK'bsel allowc*d by the 
data is 5 C The amount of crystallizatutn at three-ft»ld 
uturatum is negligible. Graphurally it can be seen tlut 
75U.f5 has a comp«)sition whx'h is close to oliviiw satu- 
ration. More olivine-normative compositions whKh crys- 
tallize un(kr equilibrium conditions are beginning the reac- 
tion olivine + lujuid -> pyroxeiw at this temperature and 
composituvn of residual liquid l lwrefore 75o V^ is in reac- 
tion relaturn with olivine at its lujuidus. and olivine is not 
observed experinwnUilly or petrographKally in the crystal- 
lization sequeiKx* It is clear tlvit 7.^lt5 1 % the crysiallizutuvn 
•'product of the celebrated 'eotcx'tie' lujutd (more properly 
' a peritcxtK lujuid sukc oIiviik is in reacturn relaturnshipk 
F(H the reasons stated in the text, we feel that this IkjukJ 


Orittin <tf litdmierout lunar hataiu 


IS 


IS derived Irom fracnonal cryslalli/alion ttf tlie parcnlal 
liquids to the Iom-K suite 

We have also given a cursory petrographic descriptHin 
of 71569 in LoMiHl ef al. (I974>. Sample 71569 tecturally 
rcsemhies 70215, one of the starting nviteriaN used to 
determine the saturatKin curves Holh riK'ks appear to have 
pyroxene before plagioclase in the petrographic cryxtalli/a* 
lion sequence However. 70215 has a bulk composituin and 
experimental crystalli/alion sequence in which these two 
minerals crystalli/e in the opposite order The petrographic 
'anomaly* was interpreted in terms of delayed nucicalum 
and irrowth of plagKKbse as a result of superuMiling dur- 
ing the crysialli/ation of 70215 On the other hand. 71569 
has a bulk composition and experimental crystallization 
sequence compatible with a petrographic sequence in 
which pyroxene crystallizes first. Noting the similarity of 
texture to 70215 it is probably that 71569 also crystallized 


under strongly superuxiled conditHins even though this 
interpretation is not required by any discrepancy hetweeu 
experimental and petrographic cryxtallizatkm sequences 
The bulk u>mpr>vitK>n of 71569 is intermediate between 
the low and high-K suites with respect to nuior elements 
Its compositKvn cannot be parental to, tw be derived from, 
the parental liquids of either the high or low-K. suite by 
fractional crystallization at the lunar surface 
Sample 71569 is also intermc*diate in crystallization 
sequence between the high and low-k suite We have stud- 
ied a high-k Apollo II basalt l(N)72 Krom petrographic 
observations, these intersertal rocks crystallize pyroxene 
before plagioclase IJamis and Ja< ksos. l97o) f xperi men- 
tal crystallization results on 10072 arc consistent with this 
observatKin The compositional and thermal interval 
between pyroxeiK* and plagKK'bse entry is greater (or 
10072(20 O than for 7l569(<iU C'k Microprobe analysis 


Table 2. Titaniferous basalt cxperinKnts 
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• All experiments performed m a mixture of Hj and C'O^ (log • 025( at I atm total 

pressure which had a pOj slightl> more reducing than the he heO equilibrium Estimated pO- 
raqged from It) " to It) '* ** over the temperature range of these experiments. .. . 

" Equilibrium pbg entry between II.V)and II.Vi. Note glassy starting material ibr (Xi-.V» (ailed 
to nucleate plag during heating or even after } days at temperature. 

f^iatural orange soil ground to uniform powder used as starting material unless specitied 
otherw ise. 

Note: Osp -= chrome spinel, qpx > quench crystals, mostly pyroxene airma! = arnulctdite. 
hiC'apx. ItX'apx » high and low calcium pyroxene. Subscripts to h'a,,, Fpb,,. and llm„ are 
(Ke Fe •¥ Mg) x UN) compositions lor olivine. arnt^<lite. and ilmcnite 
A. C. T arc wt"„ AI;Oj. OjO,. and TiO^ numbers. 
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of residual liquids 'rom these experiments conform to the 
saturation curves constructed in the earlier companion 
paper (Lostiill ef al.. 1^74). 

Orange .«*// 74220 

l-xperimental work up to 2.' kbar has been performed 
on a sample of orange miiI 74220..^? The cx>mposition of 


this starting materul in Table 1 differs slightly from that 
of 74220..' ll.SPhT. I97.t) ’which was discussed in the text 
and shimn in higs. I 4 Sample 74220.' would be capable 
of yielding the low-K ma)<K element eotcx'lic chemistry by 
equilibrium crystalli/alion of olivine On the other hand, 
our sample 74220..'7 cryskilli/es pyroxene slightly before 
plagKKlase in mild vacuum Olivine crystalli/ation yields 
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reiiidutfl liquids mmilar the strongly fraciHinulcd repre- 
sentatives of the lo«-K suite, or liquids transitional to the 
high-K suite, hut not like the parental loss-K uimpositiont. 

Small compositional dispersuins in the neighKnhiHxl of 
74220 may possibly produce parents to both high and low- 
K type major element chemistry. (It »hould he noted that 
this dispersion cannot he introduced by i>livinc fractiona- 
tion. J There is evidence that this compositHin existed at 
the surface of the Moon as a liquid. (It should, however, 
be noted that the liquid was not superheated since olivine 
phenocrysts are present.) This gbss compositKin has a wide 
lunar distribution. 

Given these circumstances, one cannot ignore the possi- 
bility that one is dealing with a primar> liquid composi- 
tion. This possibility has been recognized by Gkii n t't al 
(I974|k who have reported crystallization experiments on 
synthetic orange soil. We confirm the results of those ex- 
periments here and agree that the orange soil makes a 
better point of departure for chemical and mineralogical 
exploration of the deep lunar interior tlian d«H.‘s the Apollo 
1 1 'model basalt' or other representatives of the titani- 
(erous basalts. Assuming the primary nature of 74220. the 
experiments show olivine and orthopyroxene simul- 
taneously on the liquidus at 20 kbar. This corresponds 
to a depth to the supposed olivine pyroxene source region 
in excess of 4(M) km. Figure b indicates the migration of 
the olivine pyroxene cosaturation boundary across this 
composition space as a function of pressure. It is the 
strongly olivine normative character of 74220 which dic- 
tates such a deep sourc'e region 

Increasing degrees of melting of the olivine pyroxene 
source region at depth will produce compositional disper- 
sion of the liquids Subsequent low-pressure fractional 
crystallization of olivine in the different batches could pro- 
duce liquids which crystallize pyroxene and plagioclase in 


either order In general, the batch produced by the greater 
degrex's of nwlting al high pressure will liave pyroxene 
crystallize before plagKiclase at low pressure This property 
of crystallizing pyroxene before plagux'lase characterizes 
the high-K suite ami so if it u>mes from depth it might 
be expcx'ted to be the result of a greater degrex of partui 
nwlting than the low-K suite on the basis of its majiK 
element chemistry The trace element chemistry suggests 
the relationship is the reverse This discrepancy can be 
removed by dropping the assumption that tlw two suites 
originate at the same depth The model discussed in tlw 
text derives the high-k suite from a shallower source with 
a smaller degrex of partial irwlting The res«»lution of this 
problem dtws luit dicUite that the different depths both 
be at relatively modest depths, as was proposed in tlw 
text on otiwr grounds. 

If 74220 IS a suitable parent for oiw suite, the parent 
to the other suite is not obvious It clearly caniKit be par- 
ental to bi>th suites, if 74220 can produce titaniferous 
basalt by fractKinal crystallizatKin. it is not clear why there 
are no intermediate differentiates This contrasts with the 
continuous compositKinal variation produced by fractional 
crystallization within the low-k suite of titaniferous 
basalts. Furthermore, such a proposition is contradwted 
by the Kb. Sr. Zn. and G anomalies noted by LSPI T 
(l‘f7.H 

Alkalii'i* 

We have investigated the effect of alkali addition on the 
location of the xituration curves with results discus.sed in 
the text. Sample 70215 was used as a starting material 
to which Na^O was added as NaH( Oj. i xperiments were 
conducted at I atm in Mo capsules in a gas mixing furnace 
(Fi 2 , C'0<) to regulate pQj. Results of mwroprobe analyses 
of the residual glass arc in Fig. 6 and Table I. 
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, Differentiation of a very thick mauma body and implications for 
the source reKions of mare basalts 

David Wai.ki.r, John I.ongmi. and Javus Fmm> Hays 

Center far Earth and Planetary PhyMCA, Ht>(Tman l.ahtiT4lory. Harvard llniver\ity, Cambridge. 

Mass 02ltK 

Ab«lra<i — Mawhalance ^akulatHtns indicate that the molten layer originally covering the may 

have been several hundred kilometers deep The sol»ditic»tion of this m.*|tinn tK'ean involved some 
unusual effects resultinic from the lat ye -pressure difTercniiai in the ihiek magma layer The sunken 
cumulates from this differentiation wtmld he more iron ruh. less rcfractorv. und would contain more 
inc'ompatihle elements than wmild he expected in simple is4»Kiric fraclMm.il cryslalli/atMtn and crystal 
sinking models of magma differentiation Ihcsr lower cumulates can kive the mineialogical and 
chemical properties of the stiurce regions of the low Ii mare basalts 


Ini Kom'f' i u>n 

Skvekai. uni:s of geochemical and geophysical evidence suggest that wholesale 
melting of the outer portions of the moon occurred very early in its history, 
perhaps due to accrefional heating. The basaltic samples returned from the lunar 
missions may be partial melts of source regions which were differentiated during 
the solidification of that large magma rVean. It is the piirpi»sc ol this paper to try to 
estimate the si/e i>f this melted /one, to examine its mode of dilTeientiation. and to 
see if the pt»stiilated differentiation priulucts can have the properties of the stnircc 
regions of lunar bas.ilis. 


Mac.ma mody si/.i 

r: • f 

Kvidence for a global melting event early in lunar history will not be reviewed 
in detail here but has been deduced from petrologic, trace-element, geophysical, 
and isi>tt>pic consideralitms, l!arly identification of this event was made by Winni 
rf til. (1970) and Smith tt ol. (1970) on the basis of the complementary nature of 
mare basalts and anorthositic particles found in the soil at the Apollo II' landing 
site. The feldspathic pui tides were presumed to be characteristic of lunar crustal 
material exposed in the hmhiands. The complementary nature of the basalts and 
crustal materials was further indicated by respective negative and positive 
europium anomalies documented by many including Fhilpolts and Schnet/lcr 
( 1970) and Wakita and Schmitt ( 1970), Seismic, gravity, and rheological considera- 
tions indicated that the lunar crust w.is not of trivial thickness and must have been 
produced early lt» support mascims. Rb-Sr systematics of lunai rocks and s«»ils 
and widespread occurrence t>f ancient Rb-.Sr model ages suggested that the event 
responsible for this stale of affairs occurred very near to the formation of the 
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rm>on at —4.60^^. (PapanastassiiHJ and Wa^scrburg, 1971). These early dcduc* 
tions have been reinforced by later studies. 

The si/e and initial composition of the mugmu body are open questions at 
present, hstimates of the initial composition of the magma body or, alternatively. f 

assuming homogeneous accretion, the bulk composition of the miMin have been 
made by several authtirs; (I) Anderson <1973 >. (2) Cianuputhy and Anders <1974). 

(3) KingwtHHl and l.ssoncc ( 1970). <4) I aylttr and Jakes < 1974). .md <.^) Wiinkc e\ ul 
(1974) None of these suggested c*»mpositions appears to be capable of satisfying 
.ill geochemical, petrologic, anil gcophvsii.d constraints but sonn* models appe.ir 
to be much better candidates ih.m others. I hcie is no obvious way to make model 
(I) consistent with experimental petrology constraints on the si»urce regions of 
mare b.isalts Model (.^1 is tiH> low in silica ti* allow much pyroxene in the source 
regions of mare basalts MinJel (3) may satisfy some gross geophysical constraints 
on the pioperties of the moon, but it would be ditlicult to produce the lunar crust 
by ditferentiation of this composition. f)n the other hand, models (2) and (4) 
appear to be able to generate rocks of the lunar crust by fractional crystallization 
and give a plausible constitution for the lunar interior. Hinlges and Kushiro (1974) 
have investigated the phase relations of miniel (2) as a function of temperature 
and pressure and reached the conclusion that it could give a generally satisfactory 
account of lunar petrogenesis and lunar geophysical properties w ith the following 
exceptions. The F e/Mg ratio of lunar crustal rircks was higher than might be 
expected from model (2i and fractional crystallization of model (2) did not 
produce trirctolite as a rock type. Both these diflicultics can be overcome by 
recalcul.iting slightly more of the total I e in model (2> as oxidized I'e rather than 
metal If model (2) is oxidized it is very similar in rnajot -element composition ti> 
model <4 1 except for its lower silica content. This mtulilied moilel would 
piecipit.iie spinel dining fractional crystallization at low pressure, ie.idiiig to (’*• 
depletion of the lesidual liquid. .Since as much as of the crysi.illi/ation product 
can be spinel without piodiicing marked (’r depletion (Walker vl ill.. l97S,i.b) .md 
since the oxisteiKc of spinel lri*ctolites requires th.it there h.ive been some spinel 
.iccumiilalion. ,i successful model will precipit.ite only a minor amount oi spinel 
during dilferenliation Minor oxidation of model (2) or a slight decrease in silica in 
model (4) would be an .icceptabic solution. Therefore we have chosen hulk 
^ compositions for the large/magma system (or for the moon assuming homogene- 

ous .iccrelion) intermediate between the model of (ianapathy and Anders (1974) 
and Taylor and Jakes ( 1974). (.SiO-, 43: AM).. K.2: F'eO. 10..^; MgO. 3 1 ; CaO. 7 0.) 

These magma bulk compositions are shown on Tig. 1. Oxygen units are used to 
faciiM.iie viiliime estimates Also shown is an estimated average cimiposition for 
the liin.ii crust deduced from a compilation of published surveys of lun.ir 
feldsp.ithic rocks (Walker ef «/.. 1973) This composition is roughly equivalent to 
“anoi ihositic norite.’* C'learlv for the case wheie the anorthositic-norite crust is a 
ililleienti.ite of the huge magm.i there iiiiist be a consideiable volume of 
complement. iry material stored somewhere in the system to preset ve m.iss 
b.il.iiKi The central pioblem is to deduce the character and volume of this 
luissme iiialeii.il 
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f ig I hslim<iles of lunar hulk compt>\ition arc indicated by open cymbolt. A ic from 
Anderson <IV7.1t, (iA is from (ianapathy and Anders (IV74) calculaled »ith all Fe 
oxidi/ed: HK is from Hodges and Kushiro ( 19741 and is based on (i A with only purl of 
the Fe oxtdi/ed. RF. is from KingMxHtd and Fssene M970>; TJ is fri>m Taykn and Jake< 

<|974|. W IS from W.inke el al. i|974l fomp«*sitions v^hich satisfy petrological and 
chemical constraints on the lunar hulk composition should lie nearlv on a line hcisecen 
Mf SiO. (Mf MgO •• l eOl and the periteclic insolving spinel. »»livine. anmlhiie, and 
tiuiiKl I Walker el al 1**7 .1. |97<a.hl I he successful model comp4*sitii*ns are laK’led 
M \(iMA I he compositions labeled ('Kl'ST corresp«*nd to anorthositu norite M M l( 

Cl 'Mlll.A'I I indicates the n.itiire of the ni.iteiial coinplement.irv to the ciiist (oime«l in 
ilillcientiation of M A(iM A The diagiam is pKitted in oxygen units to fas ilitute estimates 
of the relative volumes of CKl'M and MAI If ( l^MI I All as nulls ated «*n Ihe lines 
connecting them since oxxgen compiises most of the volume of the inalerials Oxvgen 
units .ire calculated b> weighting Ihe compoiirnl lormiil.ie by Ihe number of oxygens 

they cont.iin 

Also shttwn in Pig. 1 is u set of low-pressure liquidus phase relations which can 
serve as a rough guide to crystalli/atittn behavior (Walker et al.. 1973) although the 
excess high-calcium pyroxene component (Ca: Al ’>0.5) is not considered. Crys- 
tallization al the upper surface of the magma where heat loss occurs shmild he 
ci»ntrt»lled by these liquidtis boundaries. Olivine should crystallize first (sttmewhat 
above I4IMVO followed by ci>-piecipit.ttion of spinel until plagittclase crystallizes 
at I275 (‘. Olivine and plagitKlase precipitate until I22^'(’ when low -calcium 
pyrttxene joins the crystallization. I he lunar feldspathic crust is presumably the 
non-sinking fractitm of crystalline material so priuluccd. Differential plagit)clase 
accumulation is responsible ft>r the relative buoyancy t»f the crust us well as the 
majtn comptment t>f chemical variation in the crustal n>cks. T«» explain I'c/l e ♦ 
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Mg variation in the crustal rocks we need only note that plugioclase fractionation 
t is accompanied by olivine and pyroxene fractionation which can change Fc/Fe 
Mg of the evolving liquids. Clearly this ferromagnesian material is not all 
preserved In the crust and must exist as a complementary cumulate at depth 
beneath the crust. The volume of this complementary cumulate can he estimated 
frt)m the bulk composition of the system and the seismically determined volume 
of the anorthositic-norite crust. 

Fig I illustrates how this estimate can be made. It is assumed here that the 
mafic cumulate is olivine and pyroxene, that it is devoid plaguiclase, anti that it 
has no mt>'c than y/t spinel It is also assumed th.it the bulk composition Is 
differentiated »>nly into feldsp.ithic crust and complementary mafic cumulate. The 
residual liquids t>f the system are lumped with the mafic cumulate for the present 
I'he ctunplementaiy cumulate must then represent at least i t>r 4 of the material of 
the whole magma system estimated from the volume pn*portu>ns shown m big. I 
If the complementary mafic cumulate cont.iins plagioclase or »>ther aluminous 
phases then it will represent even more of the total different!. itc Likewise, since 
there appears to be little excess high-calcium pyroxene in the anorthositic-noritc 
bulk cornpositii>n for the crust, the excess high-calciunt pyroxene in the bulk 
ci'mpositiim of the magma system must be incorporated in the complementary 
mafic cumulate. Ihis consideration would also tend tt) increase the estimated 
nniportion of materi.d represented in the complementary cumulate. .Seismic 
experiments have been interpreted tr> shtrw that the lunar crust is about f*0 km 
thick on the moon's nearside ('lt>ksb/ et al., L>74). I he ami>unt of sh.ittered 
feltlspathic crustal material versus mare basalt in the upper 20 km is uncertain and 
undmibteilly highly variable. I’or r>0 km of .morlhositic ru>rite the wh»»le system 
must extend ti> a depth of 2<Mt-27^ km depending on whether the crust is \ t»i \ of 
the system. I he 20<l-km figure agrees with the estimate I'f H»ulges and Kiishin* 
(l‘>74) who also n«>ted th.it this must be a minimum thickness for the wh»»le 
system. 

Woi'd (l^7t) h.is noted that m.iny aspects of lunar .isymmctry may be c.iused 
by a greater thickness of crust on the imron's farside If we ch.inge the crustal 
thickness estimafe to HMtkm. the whole system ctmld be us much as ^ttOkm 
thick. The presence of plagit>clase or other aluminous phases (other than the V'i 
spinel) and the inclusi»>n of the excess high-calcium pyroxene below the crust will 
further increase this estimate 

This estimate involves melting of roughly half the volume t>f the mi»on. If 
accretion t»f cold material is the heating agency to cause the melting, then the 
magma ocean cannot be very much thicker than this estimate The total 
gravitational energy released in accreting (he moon is 4(K)cal/g. which is about ^ 
suflicient (»> heat and melt the whole moon, but the accreg^im process will iclease 
heat in a manner dependent on the square of the radius f»f the grrming planet, so 
the 4(K>cal/g is not distributed uniformly . Potential energy ctrnvcrted to heat, even 
with cmnplete efficiency, will be sufficient to melt (and later superheat t subse- 
quent inf.dl only when the moon has reached i of its present radius Appreci.tble 
kinetic energy t»f accreting particles or accreli»>n of hr>t particles couUI lUsrease 
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the depth of the melted zone. Rapid accretion is necessary to prevent escape of 
the released heat. Mizutani et al. (1972) have considered the process of accretion 
and its implication for lunar thermal history in some detail. They conclude that a 
satisfactory account of the moon's thermal properties can he given by an 
accretion model which assembles cold material in about 10* yr. The depth of the 
molten layer is about 400 km, which is compatible with possible estimates of the 
macma thickness based on chemical and mass-balance arguments above. We must 
n»>w inquire how such a thick magma will solidify and what will be the character 
of the differentiation products. 


What Roc ks ari at 200-300 ksi in rm. Moon? 

Studies of the melting behavior of low-Ti. olivine-rich, mare basalts suggest 
that they would be in equilibrium with an olivine and low-calcium pyroxene 
residue at a depth of 200-3(M)km in the moon (Green et al.. 1971: Grove et al.. 
I973; Kesson. 1974; Longhi et al.. 197^; Walker et al.. 1975). We have already 
noted that the complement.ir> mahe cumulate should be composed of those 
minerals and could easily extend to this depth. An important question in lunar 
petrogenesis is whether the low-Ti marc basalt source region lies within the 
complementary mafic cumulates or in undifferentiated material. The former 
possibility has been considered unlikely by Ringwood (1975) on the grounds that 
b'e/r’e + Mg and incompatible element abundances attributed to the complemen- 
tary mafic cumulate are too low to be consistent with the properties of the low-Ti 
male basalt source region. For instance, the bulk Fe/Fe ♦ Mg of the model 
compi>sition we are considering isO. lf>; the first olivines to precipitate would have 
a F'c/I e ♦ Mg of .ibcuit t).0b whereas the source region for the low-Ti basalts must 
have a Fc/Fe ♦ Mg of 0.24-0. 29. I he crystallization of the magma at the top of the 
system where he.it loss occurs wmild precipit.ile ferromagneslan miner.ds which 
sink to ft>rm the cumulate. Hmlges and Kushiro (1974) quantitativels pursued this 
approach for a 2(M)-km magma and ci>ncludcd that about h.ilf the system would be 
composed of dunites and h.ir/burgites. These materials would indeed make 
poor-mare basalt source regions because of their low F’e/Fe ♦ Mg, low- 
inci'mpatible element abundance, and refractory nature. This is especially true for 
a perfect fractionation or t»>tally efficient differentiation model. The very si/e *>f 
the system dictates relatively slow heat loss and hence would promote efficient 
differentiation by crystal settling. However, for a magma of 4(M) km depth the very 
large size also introduces effects which reduce the efficiency of the differentiation 
process. We will now- ct'nsider these effects. 


Ei f ia IS «»i Ma(.ma SjZi: ON rll FFRHN TIATION IN THF r.FNFRAI. « AS! 

I he principal agent w hich introduces unusual effects into the diff eientialn»n of 
a veiy l.irge-'magma body is the large-'p.’v*-'‘Ure difference between top and bottom 
(-2tfkbar for 4(M)km on the moon). This pressure differential will change the 
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crystalli/ation behavior of the m;t(!ma with depth so that different minerals may 
crystalli/c. It will also increase ihe mcltinK temperature with depth at a rate such 
that the melting temperature giadient is steeper than an adiabatic temperature 
gradient. Figure 2 shows a compilation of liquidus curves of various lunar 
materials which have olivine as the low-pressure liquidus phase. Materials with 
more normative olivine have higher liquidus temperatures and the change from 
liquidus olivine to liquidus pyroxene occurs at higher pressure. The temperature 
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t ig 2. l.it|uiJiis curves as a fiinclion t»f Icmpcialure and prcs>^ure are ploUcd for lunar 
materials paving olivine on Ihc liquidus at low pressure (from t.onghi etui.. 197^, Stolper. 
1974, Walker el al , 1972. l97SM.hi. I he m.ilerials with higher Icmperalure liquidus curses 
have more normative olivine. Note lh.it there is a regular increase of the lemperaliiie of 
pyrokcne-olivinc co-saturation with increasing pressure. Olivine-pyroxene co-valuratcd 
liquids art also increasingly olivinc-normative with increasing pressure. Average satura- 
tion gradients are indicated 14072. the most noticeably deviant composition, is distinctly 
more alupiinous than the other comp«isitions. These regularities can be used to roughly 
cktrapolMte the liquidus of a composition having olivine as the low-pressure liquuius 
phase if Ihe low-pressure liquidus temperature is known 
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gradient of a magma in liquidus equilibrium with olivine would be — O.I.^^C/km; 
with pyroxene it would be — 0.4“C/km. If a magma were saturated with both 
olivine and pyroxene throughout its vertical extent the temperature gradient 
would be — 0.8"C/km. Note that this situation would be accompanied by an 
appreciable compositional gradient in the liquid since the olivine pyroxene 
liquidus boundary curve shifts to more olivinc*normativc compositions with 
increasing pressure. 

The large-magma system may be expected to crust at the top where rapid heat 
loss occurs, (irowth of the crust by heat loss through it will precipitate olivine, 
plagioclase. and pyroxene. We have noted above that much of the ferromagnesian 
material was not retaineil in the crust but s.ink due to its high density. Olivine 
which sinks from the bottom of the solid crust (at l2(Nt'(') will undergo 
compression while sinking. An adiabatic temperature rise «*f about 20'(’ would be 
experienced in sinking 4(Ht km. However, the .ictual temperature rise in the 
magma must be considerably more than this just for it to be in the liquid state 
l‘urlheimorc. as we will discuss later, the accreliorvprtHJucid thermal prtilile 
would have a ci»nsij)erable Kuly of superheated magma at depth- ('learly the 
sinking olivine will not be in thermal equilibrium with the liquid. The olivines will 
dissolve until the superheated magma is Cixded to its liquidus temperature by 
supplying the heat of fusion. Note also that this pri>cess introduced a compt>si- 
tu>nal gradient in the liquid. '1 he magma will h.ive a higher liquidus temperature 
and become mttre olivine m»rmative with depth .is olivine tiissolves. When magma 
at one ilepth is cooled to its liquidus temperatuie. olivine will no longei dissolve 
but v..in sink It* s;ij\n;ale .1 deeper level by dissolution theie. However, in linking 
thrt*ugh the saUiratetl level the coolei olivine must induce extra crysl.illi/ation t»f 
the magma as it approaches thermal equilibrium C'rystal sinking in a large magma 
may be ct*nsidered a negative im*de ol heat tr.msfer which induces cryvtalli/alittn 
at depth even thtuigh heat is lost at the lt*p of the system. 

(’ontinued crystal sinking will bring the magma at depth tt> co-saturation with 
olivine and pyroxene with a ct>rresponding compositional gradient wdh depth. 
The temperature gradient for this situation is seen in l ig. 2 ti* I'e 0.8‘(V'km. The 
ability of the system t<» reach this gravlient. rather than *‘chi*king'* and st*lidifying 
at st*me intermediate depth, ilepends on the original anu*unl of superhe.it with 
depth, the r.ite of heat U»ss at the top. and the rate of crystal settling. Rapid crystal 
addition and small superhe.it will prevent this gradient from being established 
over large ranges of depth To give a feeling of the possible magmtiule r*f the 
effect for the case where such a situatiivn is reached, we may calculate the amount 
of heat necessary \o bring an olivine sinking from the crust into equilibrium with 
the magma as it sinks t*> the bottom i>f the system. 1 his heat must be supplieil by 
crystallization of the saturated magma through which the ci'ld cryst.il sinks In 
sinking 4(M) km the olivine experiences a temperature rise t*f 2d C' due it* .idiab.ilic 
compressii>n. At 4(Mt km depth the magma is (4(MI x. h.M 20) ht*tier than the 

olivine. Using representative heat c.ipacity d.il.i fu*m Kt*bie and WaUlb.uim ( l‘W*X) 
for t*livine. abi*ut I4kcal/gfw is necessary tt* raise the t»livine to the magm.i 
temperature, b'or a heat of fusion t*f olivine »*f 22 kcal/gfw abt»ut half .is much 
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olivinc as sank into the magma would be required to crystallize In order to keep 
the olivine in thermal equilibrium with the magma. (These heat capacity values arc 
uncorrccted for pressure but serve to give a rough idea of the order of magnitude 
of the process.) This extra crystallization would, of course, not all take place at 
40(>km but would happen continuously during sinking. 

The pressure gradient also may change the nature of thcicrystalli/ing ^>hases \ 
with depth. For the case we have just discussed, the composition gradient with 
depth probably will keep olivine as the liquidus phase with depth. However, if 
there is appreciable mechanical stirring (e.g. from impacts) or convective stirring 
(to be considered later) and the compositional gradient is suppressed, the sinking 
olivine may in some cases become unstable with respect to pyroxene. .Should the 
sinking olivine armor itself with a reaction product of pyroxene, the lesser heat of 
fusion of pyroxene (l.^-IH kcal/gfw) would require that a larger number of moles 
of pyroxene precipitate than was the case of precipitation of olivine to heat the 
cold sinking olivine. The number of moles of precipitated pyroxene must be 
c»>mparable to the number of moles of sinking olivine rather than the previous 
estimate of about half as much for the simple olivine case if olivine precipitatiiui 
alone is to supply the heat At depth the reactii>n relation between olivine, 
pyroxene, and the liquid becomes a co precipitation relationship if the liquid 
comp»*sition is allowed to adjust itself according to the pressure constraints 
Therefi*re a combination of olivine and pyr.rxene precipitation would supply ihe 
necessary heat to warm the sinking olivine. In this case the sin) mg crystal would 
induce extra crystallization of the liquid by an amount from about a half t»* a 
whole of Itself, depending on the proportions of olivine and p>roxene precipi- 
tated This is to say that '« to 1 of the cry stal accumulate is material precipitated .it 
depth even though all the heat loss le.iding to this crysiallizalion occurred al the 
top of the system. 

The result of this process of heal transfer by cryst.d sinking is not only to 
accelerate crystallization at the bottom of the system but also to atfect the 
constitution of the cumulate produced. Simple efficient ditfeienti.ition mmlels not 
consuleiing this effect give sunken cumulates of diinite then harzburgiie which 
have very high Mg/Mg * l e at the bottom. The model piesenteil heie f.ivors 
harzburgites over dunites although dunites are still pivssible The cumulate will be 
less refractc'ry as a result of the introduction i>f pyroxene f-urtherinore. the 
I'e/l e * Mg of the cumulates will be higher at the bottom of the system This 
effect is illustrated schematically in Fig. The upper part of the figure indicates 
schematic liquid-solid “loops'* f«»r phases entertaining I'e-Mg solid stdutit»n \t 
lowei pressure (Pi,) the solidification curves are at lower temperiiture (7 ) than .it 
high pressure (P .«). The lower part of the figure indicates the variation t>f f e-Mg 
in cumulates produced in two different situati«ms The solid curve is fi»r fractional 
crystallization at a omsianl pressure (P..). ihe initial phase crysialli/mg .md 
sinking to the bottom to form the basal layer of the ciimiil.ite has a composition 
appropriate to the liquidus phase for the bulk composition (open symbol) .it low 
pressure. .Subsequent fractional crystallization produces extreme Iron enrich- 
ment The dashed cuive lutlicales the variation In the resulting cumiil.ite if there is 
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Fig. 3. Schematic illustration of the effect t*f pressure on the distribution of Fe and Mg in 
the cumulates produced during differentiation of a large-magma system See (ext fiK 

discussion. 


a substantial pressure yraJicnt between the lop and bottom of the melt. Crystals 
sinking at the lop of the system have the liquidiis composition as in the first case. 
However, at deeper levels the shift of the **l(H>p’* indicates that the sinking phase 
bectxmes more iron-rich as a consequence of the greater degree of crystalli/alUm. 
This is a schematic representation, since the sinking will not be isothermal, but it 
serves to illustrate that the extra crystallization induced during sinking drives up 
the I'e/I'c + Mg of the crystalline fraction of sinking maletial. I he composition of 
the basal cumulate layer will be intei mediate between the twt» vertical dashed 
lines. I he exact comptisition depends on the degree t»f l e exchange with the melt 
during sinking, the heat capacity, and heat of fusion of the phase 1 his basal layer 
will always be more iron-rich than in the firvt case. As the magma tills up with its 
settling products the pressure gradient and settling distance decrease, so the 
“loops*' at latei stages in the dilfeientiatii>n meige with the f\, **li»op,’‘ lAtreme 
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iron-eniichmcnt occurs at the end of the sequence but u smaller volume of 
iron-rich material Is produced since more iron was Incorporated in the early 
cumulate. The solid and dashed curves are shown intersecting at an Fe/Fe Mg 
value corresponding to the right vertical dashed line. This situation is possible but 
need not be the case in general. The compi'sitlon and depth of this intersection Is 
not uniquely specified by the position t»f the “loop.’* 

It should be noted this inwiease in l‘e/l e •* Mg of the basal cumulate occurs 
whether the material transfer is by convection (adiabatic) or crystal sinking (not 
adiabatic), lather way the incre.ised degree of crystallization induced during 
sinking increases the I c/l e • Mg of the crystalline material deposited at (he 
bottom of the system, ’The I'eO component of the resuliial liquid is more 
iinifoimly disti ibiitCil in the system than iliinng etlicieni ilitfercnti.ition. 

It is doubtful that the mode of solidiHcatnm discussed abt»ve would ever 
tiperate fully lor reasons to be discussed below. However, the ilisciissi«m does 
show that simple crystal settling models »»f «litfeienti.ition are quite inadequate to 
desciibe a 400-kni magm.i bt*dy .mil m.iy give incorrect insights inti» the 
constitution of the ililTeientiation products A further factor to be considered is 
that the incie.ise in (lie liqiiidus temper.iture of the system with pressure is 
considerably gie.iter (factor of 1 ih.iii the adiabatic temperatiiie gr.ulient (.iboiit 
().0^'( 'km in the lunar c.ise) I he adiabatic gradient should approximate the 
tempci.iture ilistribution in a convecting liquid If out magma system convects 
then It should solidify chietlv from the bottom up since the ailiab.it will intersect 
the steeper liquidus .it the bottom of the system first (big 4) I his mode of 
soliditicaiion will also subvert the efficiency of differenti.ition since the opportun- 
itv foi u V st. Illiquid sepal. ition by crystal settling is diminished when the crystals 
foim .d »he bottom of the system I he basal members of the solidification pioduct 
would piob.ibly be har/burgite oi Iher/olile (with garnet or spinel). A complic.i- 
tion to be noted is th.it the moifel moon composition has an appreci.ible melting 
Intel v.il st» th.it crystals max piecipil.itc over quite a dist.ince in the column .md 
settle befoie the bottom is solidified fully. Another complication, to be developed 
latet . involves likelv miti.il temper. itiiie distributions which may be gt.ivitation.illy 
st.ible ovei large parts of the m.igm.i However, neglecting these speci.il iiuti.il 
situations, the m.igma. when coo eil .it the (op. becomes exceedingly ’instable with 
respect to convection. In the presence of a superadiab.itic temperature gradient 
convection will occur in any m.tgm.i body deeper than about 10 m for rcason.ible 
values of viscosity, thermal exp.insion. and ditfusivity. The Kayleigh number for a 
400-km ni.igma with an adiabatic gradient is about I0“. The si/c outweighs all the 
other f.ictors It is clear then th.it convection should be .in import.int part of any 
analysis of magma differentiation on this scale. If solidification from the bottom 
upwards were the only process opeiating. however, the existence of .i thick lunar 
crust containing material of low l e/1 e • Mg would be iinaccounteil for 

It sluxild be noted th.it once .i material becomes partly cryst.illine it is r.ither 
effectively “bufTered” against normal convective cycling on an .idi.ib.it Should a 
parcel of m.igma with suspended crystals begin to rise, melting of the crystals 
begins which cools the parcel .is the heat of fusion is supplied. I he temperatiiie is 
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Fig. 4. The effect of having melting curves steeper than aJiahats is to induce cryvtalli/a- 
tion at the bottom of the system and during sinking Ihis effect has been discussed in 
some detail by Irvine Ascendin/convectivc magma cskiIs at ‘he top of the system 

and descends, intersecting the melting curve at depth In actual magma differentiation 
this effect is complicated by the presence of a melting interval with generally more than 

one crvstalli/ation curve 

* 

bulTcrccJ on the melting curve which makes the parcel ciK>ler than its surroundings 
since the melting curve is sleeper than !h“ adiabat. Parcels of magma which begin 
to sink experience extra crystalli/ation. Heat of fusion is liberated, warming the 
parcel. In general the density contrast produced by melting or solidification is 
much greater than the density contrast prinfuced by thermal expansion. Thcicforc 
sinking parcels continue to sink and rising parcels continue to rise even through 
they may be respectively warmer and colder than their sur^-oundings. This is the 
reverse of the temperature relations for normal convective cycling in u homogene- 
ous liquid which has an adiabatic temperature distribution. 

We have now considered lwt> “end-member” prtKcsscs which introduce 
unusual elfects into the solidification i>f large m.igmas. Both priKcsvcs will 
operate in the solidilicatiim t>f a magma to spread the heat loss sustained through 
the lop crust to the lower p.irts of the magma and accelerate solidilicalittn there.* 
The net result is a reduction in the efficienciency of differentiation by fractional 
crysialli/atiim and crystal settling because crystalli/ation will take place through- 
out the magma. We now consider a pt>ssible outct*mc t>f the solidilic.ition v>f the 
large lunar magma binly for plausible initial cr>nditu>ns when these effects ate 
taken into account. 
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To outline the nolidificalion priKcss of such a laruc-magma body we must 
know the phase relations of the material as well as the Initial temperature 
distribution. The phase relations as a function of pressure are deduced in the 
following way. The solidus is taken from Hinlges and Kushiro ( 1V74). The solidus 
temperatures turn oul to he rather insensitive to hulk comrH>sition since this 
solidus is similar to the one pn»pi»sed hy Kingwood and lAsence (IV70). The 
low-pressure liquidus temperature (olivine I is taken from putting the hulk HcO 
and MgO into koeder and Lmsiic’s (1^70) olivine-saturation model. The high- 
pressure liquidus and olivine pyroxene saturation curve can (hen be deduced hy 
extrapolation of the trend shown in fig 2 The result is shown in Fig. The 
principal heat source for the initial melting event is presumed to he accretum. The 
accretionary temperature profile shown in Fig. ^ is essentially that of Mi/utani tl 
al. ( IV72) hut has been adiusied for the larger melting interval in the present model 
composition. 


depth km 

2C : 400 600 800 



t*'ig Possible initiiil conditions fur the piimordi.il lunar moKmii ocean are illusiraled as 
disuissed in ihe text. Impi'rlani features iru.lude melting to a substantial depth, a large 
amount of superheated m:igma. and a graviialHinally stable zone ileeper than the 

temperature maximum. 
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Inspection of the initial conditions pststuiated in Fig. 5 revenU several points of 
interest. The material is molten to u depth of '•400 km. (Accretion of hot material 
would make the molten layer thicker.) .Siihscqucnt material infalling after the 
melting temperature is reached will be superheated. A decline in the accretion rate 
at the end of the accretional process allows ciHiling of the surface layer. \ 
maximum in temperature is induced In the system. The temperature profile it 
depths below the maximum is gravitationally stable and there would be no 
tendency for convective instabilities to develop. At depths shallower than the 
temperature maximum, the profile is exceedingly unstable with respect to 
convection indicating that this initial stale of .dfairs must be very transitory In 
actual fact as the rate of infall declines, the subsequent portion of .iccreting 
material will probably be falling into a convecting bath. This state of affairs or that 
which iK’curs very stum after the decay of the convectively unstable profile is 
illustrated by profile I in Fig. f>. A final feature to be noted in this differentiation 
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miHlel is the presence of a crust at the top of the system. This is impv)sed as a 
boundary condition even though there is some tendency for the crust to he 
disrupted by the convecting layer underneath. This boundary condition is 
geologically reasonable since large terrestrial lava lakes develop frozen crusts in 
the presence of convection. Furthermore, the presence of high-Mg/Mg ♦ Fe 
material in the present lunar crust is most easily explained in this manner (Taylor 
and Juke^. 1974). The problem of crust formation is discussed further in a later 
section. 

Material ciHiled and partially crystalli/ed .it the top of the system becomes 
dense enough to sink due both to thermal contraction and to the presence of 
suspended solids. It should be noted here that plagiiKlase * 2.7) should not sink 
and would probably hi>at in the fraction.ited liquids fr«>in which plagunrlase 
cryst.dli/es (calculated p ■ ^ () by the scheme of Hotting.i and M eill. I*<7(M A 
further argument that plagioclase prefeientially remains at the t*»p to be inci»rpo- 
rated in the crust is that if mote than 1^'^ »>f the sunken inateiial is plagiiKlase 
then the mass-balance c.ilculations presented earlier would require a b«Hl> of 
magma l.irger than the whole nuM>n he melted to extract the cMisi. I'he separation 
of plagioclase from the fractionated liquid and crystals which sink depletes f.u in 
this descending material 

As the cooled liquid and its mahe minerals begin to descend, further 
crystallization is induced since the melting curves (e.g. solidus) increase in 
temperature faster than ihe liquid adiabat. Kven if there is no heat exchange with 
the surroundings, the sinking material undergoes extra crystalli/.ition during 
sinking. Iherefore the actual temperature gradient followed by a descending 
convective b«Kly will be steeper than an adi ibat due to release of he.il of fiisitm. 
The gr.idient. however, will not exceed the cryst.illization cuive the b^nly folU»ws. 
It should be noted that there was a “fl»H»r“ for convective updiafts intOHliiced by 
the top of the gravit.itionally stable regime (m>te dotted line in big. h). This is not a 
“fliHir’' for downdrafts since they are colder than Ihe surroundings .iiul cont.iin 
suspended ciysials. Sinking material will descend to the bottom of the nulled 
/one The extra crysialli/alion induced during desi.enl will make the stdid niateii.d 
of higher F'e/) c ♦ Mg than .it the beginning ol Ihe descent. When the maieii.il 
slops descending there is no incentive for .iny residual liquid remaining in the 
batch to rise again for another cycle of convective overturn and further 
dilferentiation since Ihe material is now in the gravilionally stable regime Any 
further separation of residu.d liquid occurs only by upward physic. d displacement 
of Ihe liquid bv denser settled crystals 

I he thermal clfects of such a cycle can be seen in the Sequence of curves I 
through 4 tin f'ig. h. Rising of rnaten.il and ciH>ling at the tt>p of convective curve I 
will lower the maximum temperature of the system from I to 2. Sinking of the ctH»l 
m.iterial displaces the rem.iiiuler of the system upw.iids The •‘flt>or‘’ of convec- 
tive updrafts moves up from I to 2. To Ihe extent that the sinking maieiial 
exchanges heat with its surroundings. Ihe temperature shrpe of the gr.iv it.itionally 
stable legime in 2 is less than in I. The ceiling of the vonvestive zone is loweied as 
the crust at the top glows Heal h»ss out the cold bottom »«f the system induces 
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final Holidificution of sunken material and its trapped juices from the bottom up. 
Futthcr stages of this pn>ccss are illustrated hy curves .1 and 4 in Fig. f». 

The geiK'hemical novelties introduced in this differentiation scheme arc that 
Fe/Fc Mg of the sunken crystalline cumulate dtH's not show as wide a variation 
on a large scale as might he expected during efficient fractionation and that 
residual liquids on the bottom are not convectively recycled for further enrich- 
ments hut are trapped with the cumulates, ( t his material is of course depleted in 
Fu as a result of plagiiHrlase separation which iKciirs at the top of the system 
before this materi.il descends.) The first effect results from the extra crystalliza- 
tion introduced during sinking and oper.iies independently of the second effect. 
The second effect is a result <»f the acsretionary thermal nrofile. Inas much as this 
scheme of differentiation overo>mes most of the geochemical objections (King- 
wtH>d. 1975) to regarding the mafic ciimul.ites as the source material of low-Ti 
olivine-rich marc basalts, we may conclude that it is still possible to regard the 
cumulates as a reasonable source region on geochemical grounds, 


Disfissios 

The solidification process outlined in the previous section indicates that the 
bulk of lower cumulate layers produced in a large-magma b»>dy are relatively 
undifferentiated. If the priKcss operates fully the lower cumulates will differ from 
the bulk composition of the system only in having litst plagiochoe t.md l\u. etc.) to 
the crust. I he Fe/l e ♦ Mg will be roughly the same as the bulk composition and 
the full complement of Hu-depletcd residual liquid is successively inc»*rporated in 
acciimiil.ttmg layers I'his piodiict contrasts with the very magnesian dunite which 
wouM be expected in a simple fr.ictional crystallization and cry«>tal sinking model 
where hidf the volume of the system crystallizes as magnesi.m olivine and sinks 
bef»»re plagioslase precipitation initiates crust formation. Instead of a b.irren. 
magnesian dunite. a harzbiirgite forms with some layers enriched in Ku-depleted 
residual liquid component providing a more attractive source region for low-li 
mare basalts during later remelting. 

A more quantitative treatment of the process di>es not seenv possible .it this 
time. Details of the initial temperature distribution, the rates of crust formation 
and heat loss at the surface, the flow rates and patterns in the magma, and the 
degree to which chemical equilibration is achieved remain unknown. FMirthcr 
complications which requiie exact treatment for prt'per uncerstanding are the 
effects of compositional grailients introduced by piecipitating and dissolving 
material at intermediate depths. I wo features of the mtHlcl. however, seem to 
require fuither discussion first. hi*w .tie the tiace-element eniiched residu.il 
liquids priHliiced and where are they to be fiumd at the end of the solidific.ition 
process? Second, what is the process of crust loimution atop a ctmvccting. 
superheated, magma body’’ 

There is impoitant circumstantial evidence that the trapping of resi«lual liquid 
in the cumulate h.is not been complete and unii^trm. fitaniferoUs lun.ii b.isalts 
which have strong enrichments in incompatible elements (and Ti) appear to be 
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derived from source regions in the lunar interior at depths of -IOO-l.M)km. 
perhaps from a series of ilmenite cumulates derived from late residual liquids 
(Walker et al., I973a.b). If this is the correct explanation for the source region of 
high-Ti basalts then there may have been some recycling of residual liquid during 
differentiation to achieve the observed trace-element enrichments. In the 
framework of the proposed differentiation model this could happen in several 
ways Crystal sinking may to a certain extent physically displace residual liquid 
upward from the lower cumulates. Alternatively the solidification by heat loss 
thrmigh the bottom can have an effect analogous to /one reOning As the solidified 
zone advances upward through the lower cumulates, the incompatible elements 
may be driven upward if there is good intercumulus plumbing. Kither w.iy the 
Fe/Kc ♦ Mg t»f the residual liquid is ’“buffered** by the cryst.illine phases of the 
cumulate while the degree of incompatible element enrichment increases. Finally 
it should be noted that as substantial degrees of solidification of the system arc 
obtained, the pressure differential in the smaller voiume of residual magma is less, 
so that the effects described here beci*me less imptirtant. I<» the extent that the 
unusual size effects do not totally contnd the differentiation, the remaining magma 
fraction will differ from the original bulk ct»mposition and later cumulates mav 
increase in F'e/F'c - Mg .Specifically, curve 4 |n Fig h probably does not 
correspond to a real situation because changes in the liquid will have been 
introduced which modify the phase relations and temperature profiles. The thick 
magma effects will be mi»st effective in the early stages of differentiation. 

The formation of a crust on a large Kniy of superheated magma also requires 
further inspection FUistcnce of a crust was assumed as a boundary camdition in 
the discussion ab»>ve because the mtH»n has a crust and its relatively high«.Mg/re 
ratio IS most easily explained if it is produced by freezing the top of the magma 
rather than by later differentiation from mateMal beneath F'or crustal formatnm to 
occur in this way it is necessary that ascemiing magma p.irccis lose enough heat at 
the top of the system t» crystallize at least to the point of ptagioclase pres ipitatum 
If the feldspar is fiec ti* physically separate from the c»»oled .nul paitly 
crystallized material which then sinks, this feldspathic segregate becomes the 
crust To the extent that feldspar tun separate, the crustal material will be 
gravitationally stable against sinking Ip for anoithositic norite 2.9 whereas p ft»r 
the liquid ^ 3.0) It remains to be sh»>wn that enough ciMiling of the superheated 
magma can occur at the top tif the convective zone to .illow partial cry stallization 
and feldspar separation before the magma has substantially solidified from the 
bottom upward 

Details of the initiation of crust formation are not clear t>ut may resemble the 
following. Initially the heat flux ciHiling the tup of the superheated magma is 
radiaiKin into space and has a magnitude 7 cal/cm^ * sec for assumed “bl.ick 
biHly'* behavior This flux is sufficient to ciHil and crystallize the magma al the rate 
of lOm/day This rale i»f crustal growth cannot be maintained after solidification 
of the top layer limits the heat transfer to the conductive miHje. Furthermore this 
rate t>f crustal growth is uhv rapid tu allow concentration of plaguK'lase in the crust 
by differential buoyancy. I hercfore the rapidly frozen crust will sink if disrupted. 
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as it most certainly would he by continued accrctionary bombardment from above 
and turbulence in the maitma below. Crustal fragments would tend to dissolve 
rather quickly duriig sinking into the superheated magma beneath. .Since the crust 
at this stage is undilTerentiated “chill margin.” this dissolution does not introduce 
compositional gradients but does serve to reduce the temperature of the magma at 
shallow depths beneath the rapidly reforming crust. As the rate t»f crust growth 
decreases during the change to ccmductive (rather than radiative! heal loss at the 
top and as the thermal gradient decreases at shallow depths in the magma, it 
becomes possible to separate plagioclase and olivine. Upw.ud plagioclase con* 
ceniration leads to the establishment of .1 crust which no longer sinks when 
disrupted. This insures (hat the slow conductive heat loss, condui/^ve to plagiiK*- 
lase separation and further crust.il growth, continues. The important feature r>f 
this proposed mode of crust formation is that the ci>oling power of radiative heat 
loss cun be distributed to shallow depths by crustal resorption A relatively 
sl.)w -cooling regime which is then established directly beneath the “chill margin** 
allows olivinc-plagioclase separation and formalii»n of a non-sinking crust. 1 his 
non-sinking crust is then stable against perturbations and grows by ditferential 
plagioclase retention from magma rrystalli/ing slowly at its base by conductive 
heat loss. Kuptures in this crust shmild heal by a prtKess like that leading to 
crustal formation when hot magma i> again exposed to radiative ctHtling at the 
surface. 

For the proposed process to be successful, the stable crust must be established 
before large quantities of the partially ciMtled layer beneath the crust are 
subducted by convection. It is possible that the proposed process may succeed. 
Crystal separation rates leading to crust formation (subtraction of olivine by 
settling) are driven by liquid-solid density contrasts and may be more rapid than 
the ctmveclion rates ilrivcn by therm.il contrasts in the liquid alone which act to 
reduce the opportunity for establishing a /one conductive to plagioclase-olivine 
separation. 

At knoHlt di/nuntx — Wc h.«vc h«;nc(ilcJ frttm JiskUsMons with S. K. Kesst^n. M J. O'Mai . 1 . .mJ A K. 
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here Wc Ih.tnk N Jaiks4in for assisianve in manuscript prcpurutuMi rbis work was supp«>ricJ by 
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ABSTRACT 


Experimental crystallisation of a lunar plcrlte composition 
(12002) at controlled linear cooling rates produces systematic 
changes In the temperature of appearance of the crystalline phases, 
the texture, and crystal morphology as a function of cooling rate. 
Phases crystallise In the order olivine, chrooe^splnel, pyroxene, 
plagloclase, llmenlte during equilibrium crystallization but 
llmenlte and plagloclase reverse their order of appearance and 
silica crystallizes In the groundmass during controlled cooling 
experiments. The partition of Iron and magnesium between olivine 
and liquid (K^ ■ 0.33) Is Independent of cooling rate (0.5>2000*C/ 
hour), temperatuz*e (1325-600*0, and pressure (0-12 kb). Compari- 
son of the olivine nucleatlon densities In the lunar sample and In 
the experiments indicates that the sample began cooling at about 
l*C/hour. Pyroxene size, chemistry, and growth Instability 
spaclngs ("swallowtails”) as well as groundmass coarseness all 
suggest the cooling rate subsequently decreased by as much as a 
factor of 10 or more. The porphyrltlc texture of this sample, 
then. Is produced at a decreasing , rather than a dlscontlnuously 
increasing, cooling rate. 




INTRODUCTION 

Plcrltle basalts have been sanpled by the Apollo 12 and 15 
lunar Missions ; these rocks show a variety of textures and grain 
sizes ranging from vltrophyrlc to porphyrltlc to gabbrole. Com- 
parison of the compositions of experimental llquldus phases with 
the most Mgneslan core compositions of the natural examples of 
those minerals have shown that some of these porphyries (Including 
those with vltrophyrlc groundmass) crystallized from a magma 
batch of composition equivalent to the hand specimen composition 
(Green and others, 1971; Grove and others, 1973). 

Much controversy has surrounded the possible Interpretations 
of magma cooling history deduced from the two-stage crystalliza- 
tion history (phenocrysts then groundmass) Implicit In porphyrltlc 
textures. Is a two-stage cooling history the proper Interpretation 
of a two-stage crystallization history? It Is the purpose of 
this paper to reconstruct quantitatively the cooling history of 
12002, a porphyrltlc plcrlte. In an attempt to answer this 
question. 

EXPERir<£NTAL METHODS 

In our previous studies of 12002 (Grove and others, 1973) > 
molybdenum and Ingot-grade Iron capsules were used for equilibrium 
experimental crystallization. These materials interacted with the 
silicate charges to a certain extent resulting In significant Iron 


loss. It has since been shown (Kesson, 197Aa,b; Walker and others, 
1975) that Iron loss can be virtually eliminated by use of very 
high purity Irun (<20 ppm total impurities) for lunar high- and 
low-Tl mare basalts. Presumably the Iron loss experienced In the 
previous work In Ingot-grade Iron capsules In scaled silica tubes 
resulted from reduction of some iron In the charge by Impurities 
In the Iron of the capsule. Carbon and phosphorous are the most 
likely agents of this reduction (M. J. O'Hara, personal communica- 
tion). Samples equilibrated In high purity Iron containers will 
not gain or lose iron if the material Is Iron-saturated. This 
appears to be the case for many lunar basalts. Examples of lunar 
basalts which are not Iron-saturated exist (see appendix discussion 
of Walker and others, 1975) and these materials gain Iron when 
equilibrated. in iron capsules. Evidence that rock 12002 was 
Iron-saturated during Its crystallization Is given by Grove and 
others (1973) who reported iron spherules Included In 12002 
olivines. In addition we have observed Iron In pyroxenes and In 
the groundmass of 12002 as well, though it Is not abundant (<1X). 
Therefore use of iron capsules for 12002 should be successful in 
maintaining constant iron content. Conversely, If constant 
Iron content Is maintained, this confirms the Inference that 12002 
was Iron-saturated during crystallization. This appears to have 
been the ease. 

The present equilibrium crystallization experiments on 12002 
were conducted In high purity Iron. The starting material Is the 


SAM devltriried glass used by Grove and others (1973). A single 
capsule was reused within different evacuated silica tubes for 
experiments In a mild vacuum according to procedures discussed In 
Longhl and others (197A).* a fresh charge was loaded for each 
tcaperature Investigated. Capsule recovery vas made possible by 
filling only a saall part of the capsule. The evacuated void 
space apparently caused partial alkali extraction froa the cnarge. 
High pressure equilibration In Iron capsules was done according . 
to the procedures discussed In Walker and others (1972) In a 
piston-cylinder apparatus. Drying of furnace asseablles and 
charges was done at 1100*C In lOX-N^ 90S rather than pure 
for 15 minutes rather than 30 minutes. Tungsten-rhenlua (W3S-W25S) 
thermocouples were used exclusively In the high pressure experi- 
ments with Iron capsules. 

Experiments were also performed In high purity Iron In 
evacuated silica tubes at programed cooling rates to study the 

*The effects of Insufficient care to prevent capsule oxidation 
can be seen In the appendix. Table 1. Runs il88 and #18<I were the 
last experiments performed with the capsule and were done several 
months after the other runs. The acetone In the storage bottle 
had nearly evaporated. Although no obvious rusting of the cap- 
sules had occurred, alcroscoplc oxidation of the capsule Is the 
most probable explanation for the aodest Iron gain these nins 
have experienced. 


cooling history of 12002. The capsules were prepared from 5 mm 
rod stock. All machining was done under heavy oil and capsules 
were stored In oil or acetone prior to loading and use (Longhl 
and others, 1975). The O.D. was turned to A. 5 an and the I.D. 
drilled to 2.5 ma. Lids were friction fitted plugs. The cylin- 
drical cavity, which was filled with the charge, was 2.5 me In 
diameter and A ma deep. Capsules In evacuated silica glass tubes 
were suspended In a Pt-Rh wound resistance furnace. Temperature 
was monitored and controlled with Pt-PtlORh thermocouples which 
were calibrated against the melting points of gold (106A*c), 
lithium aetaslllcatc (1203*C), and dlopslde (1391*C) charges 
used In the same configuration as the run charges. Two different 
furnaces were used: a ssialler model (1* I.D.) holding a single 

charge and a larger bore model (2-1/8" I.D.) holding up to 3 
charges. Experiments which were cooled at faster than l"C/hour 
were first equilibrated at 1350 ♦ 5*C for 8 hours or more. All 
cooling experiments were started at 1350 i 5*C to give a liquid 
of uniform structure as the starting material for all experiments. 
Charges were then cooled at a specified rate to a specified 
temperature and then quenched In water. Examination of the products 
was made by polishing an axial section through the cylindrical 
charge for reflected light microscopy and electron probe analysis. 
Capsules were not reused. 

Programming of the temperature drop In these experiments 
was done two ways. A Eurotherm D120W electronic temperature 



protraa unit Installed between the theraoeouple which sensed 
the furnace teaperature and the control/power-supply package which 
regulated the furnace power in response to the theraocouple signal. 
The program unit generated a linearly varying potential which 
caused the regulation systea to lower the temperature at a rate 
corresponding to a linear change in thermocouple EN7. This does 
not correspond strictly to a linear teaperature drop for a Pt- 
PtlORh theraocouple but the deviations are negligible for the 
resolution of this study. For fast cooling rates (300-2000*C/hour) 
the program and control systea was bypassed and furnace power 
sinply shut off. Cooling proceeded at a rate determined by the 
thermal properties of the furnace. Two furnaces of different 
design and insulation were used to allow some variation of cooling 
rate In this rapid cooling regime. The cooling rates produced 
are not linear in this situation and' decrease with tine. Dotted 
lines on Figure 2 indicate the approximate cooling trajectories. 

The values quoted for cooling rate are simply the change In tem-- 
perature divided by the length of time elapsed to produce the 
texperature drop. Once again the deviations from linearity (a 
factor of 3) are inconsequential to a study of this resolution 
covering cooling rates which vary by a factor of kOOO (0.5- 
2000*C/hour) . 

Iron metal is present in all experiments as the result of 
the use of iron capsules. Microprobe analysis of equilibrium 
run products shows the Iron loss or gain is negligibly small. 

It is, however, quite clear that iron metal is not an inert 


container. The silicate charge corrodes the contsincr, estlng 
embayments into the capsule walls and bottom. The iron so dis- 
solved Is reprecipitated as prominences on other parts of the 
capsule or as iron globules in the sllieste charge. Iron apparently 
dissolves and repreclpltates even though the Iron content of the 
silicate charge remins unchanged. The same effects were seen in 
the controlled cooling runs but the presence of strongly xoned 
silicates makes It difficult to estimate the amount of iron gain 
or loss In the systea. Considerstion of llquidus olivine chemistry 
(below) suggests that net gain or loss of Iron Is also small in 
these experiments. Iron repr.eclpltated In the controlled cooling 
experiments does not seem to show much variation in morphology 
with cooling rate. Experiments cooled to the solidus do not show 
noticeably more Iron transfer than experiments quenched from near 
the llquidus. The transfer of iron nay happen mostly during the 
Initial equilibration above the llquidus at 1350*C. 

EXPEhIMENTAL RESULTS (EQUILIBRIUM ) 

The equilibrium phase relations in high purity iron are 
shown in Figure 1 as a function of temperature and pressure. 

The difference in phase stabilities from those reported in Grove 
and others (1973) are small. The principal change at low pressure 
is that pyroxene crystallisation has started at 1175*C rather than 
1225*C and the olivine llquidus temperature is about 10* higher. 
These changes are readily understood as a result of the iron loss 


In th« earlier ezperlaenta. Reduction of the total PeO^MgO 
available to sake rerroaa£neslan silicates such as olivine 
depresses the olivine llquldus, even though Mg/Pe-*^ng of the sili- 
cate Increases. Silica activity Increases, therefore pyr‘>zene 
crystallisation Is enhanced. At high pressures the principal 
change Is that spinel Is no longer present at the olivine-*’ 
orthopyrozene llquldus pressure. The Ho capsules probably sta- 
bilised spinel with a very small addition of ozldlsed nolybdenusi 
which Is partitioned ezcluslvely Into the spinel. The ollvlnev- 
orthopyrozene llquldus pressure Is 12-1/2 kb rather than 15 kb. 

* The most noticeable differences between the previous and 
present ezperlments Is to be found in the higher Pe/Fe^Hg of 
the phases In the present ezperlments. llquldus 

olivine at pressures up to 12-1/2 kb, whereas iron loss in the 
earlier ezperiments produced olivine as magnesian as POgg. At 
low pressure, pyrozene crystallization begins with Vo^Engg In high 
purity iron as opposed to earlier ezperlments. 

In general the phases present are unzoned and show only minor 
heterogeneities. 

Details of run products and phase compositions can be found 
in Tables 1 and 2 and In the appendiz. 

EXPERI^NTAL RESULTS (CONTROLLED COOLING ) 

Ezperlcental crystallization of 12002 at controlled cooling 
rates was performed In order to Investigate the variations 


observed In changing the cooling rate and the teziperature from 
which the ezperiaent was quenched. The latter parameter was 
chosen in order to arrest the crystallization process at various 
stages. Other possible variables which might affect the crystal- 
lization behavior (such as large superheating or episodic changes 
of physical conditions such as pressure) were not Investigated. 

Plgurrn 2 shows the results of these ezperiments as a function 
of these two variables. All ezperlments were started above the 
llquldus at 1350 * 5*C and cooled to the temperature Indicated 
by the circles at the cooling rate corresponding to the value of 
the abcissa, then quenched. Phases crystallize In the order 
olivine, spinel, pyrozene, llmenlte, plagioclase (and silica). 

The equilibrium crystallization sequence of 12002 under the same 
conditions is included in Figure 2 for comparison. It Is to be 
noted that the order of crystallization of olivine, spinel, and 
pyrozene is the same but plagioclase and llmenlte have reversed 
their order relative to one another compared to the equilibrium 
sequence. The curves labelled "olivine in” etc. indicate that all 
ezperiments at temperatures lower than the curve have crystallized 
olivine, etc. Figure 2 shows that there Is a systematic decrease 
in the temperature at which each phase appears with increase In 
cooling rate. The eztent of this decrease is largest for the late 
crystallizing phases; i.e., the curvature is greatest for the 
lowest curves. (The top curves for olivine and spinel may not 
be observably depressed. The ezperiswnts on which this very 
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•light depression Is based report the theraocouple tenperature 
of the Bost quickly cooled run. While the thermocouple was 
cooling at *2000*C/hour It Is conceivable that the teaperaturc 
of the Bore icasslve charge* In Its silica tube did not fall quite 
this fast.) 

The olivines, spinels, and pyroxenes crystallise as soned 
alcrophenocrysts. Changes In aorphology and chealstry are dls* 
cussed below, llaenlte crystallises as laths. In all experiaents 
where llaenlte has crystallised, Cr-splnel phenocrysts have 
developed nantles of Cr-ulvOsplnel (see Table 2 of appendix). 

This Bay be suggestive of a reaction relation between liquid and 
Cr-splnel which produced Cr-ulvosplnel and llaenlte; however, 
no other textural relations observed support this possibility, 
llaenlte Is not observed In contact with nor In overgrowth relation 
to the nantled spinels. It may be that the reaction which mantles 
the Cr spinels with ulvosplnel Is Independent of llaenlte appear- 
ance but coincidentally occurs at about the same temperature at 
each cooling rate. Nehru and others (197^) have concluded on 
the basis of slollsr relations exhibited In Apollo 15 basalts 
that pyroxene participates In this reaction. In experiments 
cooled at slower than 50*C/hour a groundmass of plagloclase, 
pyroxene, and silica crystallites. At faster cooling rates no 
groundnass crystallizes above 600*C where our observations end. 

The grain sice of the groundaass increases with decreasing 
cooling rate. 


COMPARISON WITH EQUILIBRIUM CRYS T ALLIZATION 

The crystallization Ir.'* 'r.ed during controlled cooling of 
12002 differs In three important ways from t( e equilibrium 
crystallisation. Temperature of Initial phase appearance Is 
depressed, phases produced are chemically soned, and metastable 
phases crystallize. These three differences are Interrelated. 

The fact that the temperature at which each phase appears 
decreases with Increasing cooling rate may be partly due to the 
higher probability of nucleatlon near the equilibrium tenperature 
at the slower cooling rates. However, th^s depression of phase- 
appearance tenperature Is also undoubted y aided by growth of 
soned phases, which causes substantial departures of the composi- 
tion of the residual liquid from the equilibrium composition. For 
example, fra'ctlonatlon of augite (rather than the equilibrium 
plgeonite) depletes the melt in Ca while surface equilibrium 
growth of soned olivine and pyr'>xei»t enriches the Pe/Pe*y*g of 
the residual liquids (to 0.9 v& . 0.6 at bulk equilibrium, see 
Figure 7). These departures fr^n bulk equilibrium suppress 
plagloclase and enhance ilraenlte growth potentials. Together 
with the possibility of delated plagloclase nucleatlon, they are 
responsible for the observat on that ilmenite and plagloclase 
have reversed their order of crystallization relative to the 
equilibrium sequence. It should be noted that the growth of 
metastable phases is also related to departures of the residual 
liquid from the equilibrium composition. Silica and ulvosplnel 



do not have a atablllty field betMeen the equllibrliai llquldua 
and aolidua although ulvSaplnel nay have a aubaolldus atablllty 
field. The growth of ailica in the groundaaaa occura beeauae 
the residual liquida are ailica enriched and have high Pe/Fe>Hg 
ao that equilibrium of a ailica phaae with olivine ia poaaible. 

It appeara that groaa errors may result in trying to use 
equilibrium crystallisation data to asaign tenperatujres to 
petrographically obaerved cryatallisation sequences in extrusive 
rocks or rocks from small Intrusive bodies. This can be the 
case even if the order of phase appearance experimentally matches 
the petrographically observed sequence. The error is more severe 
for more quickly cooled rocks. 

TEXTURAL VARIATIONS 

Systematic variations in crystal morphology with cooling 
rate are observed. All textures are porphyritic. At cooling 
rates in excess of 50*C/hour run products consisted of skeletal 
phenocrysts of olivine, Cr-splnel, pyroxene and ilmenite surrounded 
by glass. At greater than 500*C/hour olivines have a platy 
skeletal habit with elaborate and regular cross-ornamentation 
between the faces of the plates. Maximum dimensions of the plates 
can be of the order of 1 mm but the thickness of the plates does 
not exceed 30 um. Spinels occur as tiny specks or minute den- 
dritic crosses of the order of 5-10 um. Pyroxenes occur as 
elaborately ornamented dendritic sprays which form masses as large 


as 300 um although the individual pyroxenes rarely have a thick- 
ness greater than 30 vm. Figure 3 shows an experiment cooled at 
500*C/hour which exhibits these features. Also shown in Figure 3 
at the same scale is a photomicrograph of 12009. a vitrophyre 
with composition quite similar to 12002. The striking textural 
similarities suggest that this portion of 12009 may have cooled 
at about 500* C/hour. 

Figure A shows repi*esentative textures at slower cooling 
rates. As noted above, and shown in Figure Aa. skeletal olivine, 
spinel, and pyroxene occur in a glass matrix for experiments 
cooled at rates greater than 50*C/hour. The skeletal olivines 
are not as flattened as at 500*C/hour and the dendritic pyroxene 
is clearly coarser. Figures AB. C. and D show the products of 
progressively more slowly cooled runs. In B the phenocrysts are 
set in a microcrystalline plagioclase-pyroxene groundcass. 

Silica is pi^bably present. AC and D show progressive coarsening 
of this groundmass with silica present. AB, C. and D alao show 
progressively coarsening olivine and pyroxene phenocrysta which 
are more equant than at faster cooling rates. Some evidence of 
skeletal morphology occasionally remaina even at the slowest 
rates . 

Cr-spinel also shows systematic morphology changes with 
cooling rate as seen in Figure 5 . At all cooling rates, spinel 
occurrence strongly suggests nucleatlon on other preexisting 
solid phases as does no other phase observed. At faster than 
500*C/hour (Figure 5A) minute apinel specka are intimately distrib- 
uted within the ornamentation of skeletal olivines. At faster than 
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50 *C/hour spln«l occurs ss alnutc dendritic crosses on the ed(es 
of olivine phenocrysts (Figure 5b). At slower coollnc rates the 
Intlxiste occurrence of spinel with olivine Is replsced by the 
occurrence of spinel with Iron eetsl (Figure 5C). Equant, euhedral 
spinel Invariably grows on the capsule walls or on iron spherules 
In the charge. The aorphology and occurrence of Cr-splnels with 
ulvosplnel aantles grown at the slower rates observed aost nearly 
correspond to the spinels In 12002 (coapare Figures 5C and 0). 

It should be esphaslted that there Is a range In aorphologles 
of all phases grown at any particular cooling rate and therefore 
only a crude estimate of cooling rate can be based on crystal 
Borphology alone. For example, some large equant, faceted olivine 
phenocrysts In 12002 strongly resemble some olivine phenocrysts 
grown experimentally at a cooling rate of ^0.6*C/hour. These 
crystals arc of comparable slse and morphology so It Is tempting 
to say the olivines In 12002 grew at this rate. This Inference 
Is Insecure, however, since quite a range of olivine sites and 
morphologies are present In both 12002 and experiments at this 
cooling rate. It would, however, probably be safe to say that 
12002 cooled more slowly than 500*C/hour on the basis of olivine 
morphology since there Is no overlap of aorphologles In 12002 
with those grown so fast. 

PORPHTBITIC TEXTURE (1 ) 

We are now In a position to give a partial answer to the 
question posed in the introduction about the interpretation of 


porphyrltlc texture. These experiments, employing s single 

cooling rate, produce porphyrltlc textures. Therefore, two stage 

crystallisation Implicit In a porphyrltlc texture does not 

necessarily laply a two-stage cooling history. This was clear 

to Bowen (191^ » p. 257). Lofgren and others (197^) hsve also 

produced porphyrltlc textures at single cooling rates In quarts- 

norsiatlve, pyroxene-porphyry compositions under different experl- 

aentsl conditions for containment and p,. control. Such behavior 

”2 

may be a rather general property of basaltic compositions. This 
does not negate the possibility that some porphyries do have a 
two-stage cooling history; however. It cannot be safely Inferred 
that they do on the basis of porphyrltlc texture alone. 

A cooling history eoeaionly Inferred froa the crystalllxatlon 
history recorded In porphyrltlc texture Is that "slow* cooling 
is responsible for growth of phenocrysts, after which "fast" 
cooling Is responsible for growth of the groundmass. The two 
stages of the cooling history may be separated by a discrete 
event such as eruption of the magma. We previously Inferred 
such a cooling history for 12002 (Grove and others, 1973) on the 
basis of details of pyroxene chemistry as well as porphyrltlc 
texture. Pyroxene Is a major phase In both the phenocryst and 
groundmass population and so Its chemical variations record much 
of the crystallization hlstoz^ of 12002. This general approach 
has been extensively employed in other Investigations by many 
authors (e.g., Bence and others, 1971; Hollister and others, 1971). 
Pyroxene phenocrysts In 12002 have Tl/Al of near 1 /A which is the 
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value observed In equlllbrlun crystallization of 12002. In 
contrast pyroxenes of the groundaasa haveTl/Al greater than 1/2, 
values clearly not consistent with the observed equlllbrlun 
crystalllk.atlon. The two stages of crystallisation were then seen 
as quasl-equlllbrlun crystal growth (phenocrysts) and non- 
equlllbrlun crystal growth (groundnass) . This la just the sort 
of behavior to be expected In a two-stage cooling history such 
as we previously Inferred for 12002, If one associates quasl- 
equlllbrlun growth with "slow" cooling and non-equlllbrlun growth 
with "fast" cooling. We have noted above that this two-stage 
cooling history Is not a Mandatory deduction based on porphyrltle 
texture alone. We oust now Inquire whether the evidence fron 
pyroxene chemistry combined with porphyrltle texture Is sufficient 
to sake that deduction mandatory. 

PTFOXENE CHE:*ISTRY VARIATIONS 

Pyroxene compositional variation Is shown in Figure 6 for 
controlled cooling experiments, equilibrium experiments, and 
12002 natural pyroxenes. In both natural and controlled cooling 
experlaenta, phenocrysts zone from plgeonlte cores to auglte 
rlas with sectorial complications. In both natural and experi- 
mental phenocrysts Tl/Al Is close to 1/k; both T1 and A1 Increase 
from oore to rim, reflecting Increasing Incorporation of T1 and 
A1 with decreasing temperature, which is also seen In the 
equilibrium sequence (Table 1). Both natural and controlled 
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cooling groundnass pyroxenes are of different composition from 
their respective phenocrysts and arc sub-calcic with variable 
Iron-enrlehnent which can be extreme. Tl/Al can be greater than 
1/2 in the groundmass of both natural and controlled cooling 
crystallisation products. In contrast to the pyroxenes grown In 
controlled cooling experiments, the very limited major element 
compositional variation, the lack of Iron-enriched late pyroxenes, 
and the lack of Tl/Al > 1/2 In the equilibrium pyroxenes Indicate 
that equilibrium crystallisation Is not a good model for the 
natural pyroxene crystallisation. 

Although the controlled cooling experiments shown In Figure 8 
provide a better approximation to the natural pyroxenes, there 
are still noticeable discrepancies. The experimental pyroxenes 
show larger amounts of A1 and Tl, and, to a modest extent, Fe, 
than the natural counterparts. The experimental phenocrysts have 
maximum Tl^Al a factor of ^2 greater than the maximum seen In the 
natural pyroxenes. Experimental auglte margins are slightly Iron- 
enrlched and experimental groundmass pyroxenes are markedly Iron- 
enrlched compared to the natural counterparts. 

The highly skeletal pyroxenes grown at faster than 25*C/hour 
have the highest Tl-fAl observed. The equilibrium pyroxenes have 
the lowest Tl't'Al observed. If equilibrium crystallization occurs 
at an Infinitely slow cooling rate then the amount of Tl^Al 
Incorporated In the pyroxenes would appear to correlate crudely 
with cooling rate since pyroxenes grown at l^teraedlate rates 
have maximum Tl-fAl Intermediate between the equilibrium and most 




r«pld eoollns values. The reason for the lesser asiount of Tl^Al 
In 12002 pyroxenes is probably that pyroxene In rock 12002 gree 
at a rate slower than our controlled coollnc experiments. The 
Ti^Al enrlchxent is thought to result In part from the delay In 
growth of phases which deplete T1 and A1 from the melt. We noted 
above that the delay in iljsenite (with ulvSspinel) and plagioclase 
appearance increases with cooling rate as seen in Figure 2. 
Furthermore, during rapid cooling T1 and A1 may not have sufficient 
opportunity to diffuse away froai the growth surface and may become 
incorporated into the growing pyroxene. These ^wo processes 
produce the same qualitative result; more rapid cooling leads to 
greater Ti^Al enrichment. The relative importance of the two 
effects is not known. 

The lesser amount of iron enrichment in the natural 12002 
pyroxene may also be explained as a result of somewhat slower 
cooling rate. Fractionation of Fe/Mg by surface equilibrium 
growth of zoned phenocrysts will be less extreme at slower 
cooling rates as dlffusional exchange with phenocryst cores is 
more effective. The limit of this process is seen in the trivial 
iron»enrichaent of the pyroxenes grown under bulk equilibrium 
conditions. 

The final feature of the pyroxene compositional variation 
remaining to be accounted for is the change in Tl/Al from 1/A 
in the phenocrysts to>l/2 in the groundmass. This phenostenon 
is basically a depletion of A1 at relatively static Ti content. 

We have already noted that plagioclase entry is likely to put 
a "lid* on the A1 content of the liquid and hence of the pyroxene. 

If plagioclase is supersaturated with respect to the non-equilibrium 


t^sidual liquid a large amount of plagioclase may crystallise 
all at once upon nucleation. This would cause an essentially 
discontinuous subtraction of A1 from the melt which is reflected 
in the abrupt drop in A1 content of the pyroxenes. Since the 
groundmass pyroxenes seem to crystallize penecontenporaneously 
with the plagioclase, some local areas of the melt may not be 
depleted by prior plagioclase crystallisation whereas others are 
depleted when the pyroxene grows there. This accounts for the 
fact that some of the groundMss pyroxenes in the experimental 
ease still have Ti/Al of 1/A while some have >1/2. In the 
natural case we infer that the residual liquid was already 
uniformly depleted by plagioclase crystallization before the growth 
of groundmass pyroxene since the 1/A trend was not observed. 

PORPHYRITIC TEXTURE (II ) 

We have now considered pyroxene compositional variation 
produced in one stage cooling experiments. The production of 
phenocrysts with Ti/Al the "quasi-equilibrium" value of 1/A 
and groundmass pyroxenes with the "non-equilibrium” value of 
>1/2 has been accomplished at constant cooling rate. Clearly 
then the evidence of Ti/Al variation, even coupled with the 
evidence of porphyritic texture, is insufficient to make a two 
stage (slow then fast) cooling history a mandatory deduction 
from a two-stage crystallization history. Supersaturated crys- 
tallization, which can be deduced from the Ti/Al of the pyroxenes 
cannot be safely equated with fast or slow cooling. It appears 


that taitur* and pyroien* coapoaltlon cannot ba uaad to tranalata 
unasblguously eryatallltatlon history Into cooling history. 

OLIVINX CtgXlSTKY VARIATIONS 

Ths variation In ollvlns chanlstry Is considerably less 
than that shoen by the pyroaenes. This Is a result not only 
of the Bore Halted solid solution possible In olivine but also 
of the fact that olivine Is the llquldus phase and hence there 
Is no fractionation of solid phases prior to olivine which night 
affect the coaposltlon of the liquid at the start of olivine 
crystallization. Olivines crystallised In controlled cooling 
ezperlaents are extensively zoned at all cooling rates. Zoning 
Is not necessarily regular or concentric, especially in the 
skeletal aorphologles grown at high cooling rates. Figure 7 
shows the range of olivine compositions seen in controlled cooling 
experiments as a function of terperature. Only the Fe/Pe«^Rg 
variable Is considered since any variability In Ca, Al, T1 or Cr 
was saall or below the observable Halt of detection. Experlnents 
for all cooling rates are Included In Figure 9 along with the 
eor.posltlons of the equlllbrlun olivines as a function of tempera- 
ture. The olivines grown In controlled cooling experiments have 
their most fayalltlc olivines (rims) Halted by the equilibrium 
value for the composition of the residual liquid, regardless of 
cooling rate, at temperatures between the equilibrium llquldus 
snd solidus. At temperatures at least 100*C below the equlllbrlun 
solidus the olivines continue to maintain surface equilibrium In 


terms of Fe/Hg distribution with the residual liquid. This is 
seen In Figure 6 where the most fayalltlc olivine found In an 
experlnent Is plotted against the Fe/Fe-*^ 'S of the adjacent residual 
liquid from that experlnent. There Is very little deviation froa 
a Kq value (Roeder and EmsHe, 1970; Longhl and others, 1975) 
of 0 . 30 - 0 . 33 . The saxM values are also found for the equlllbrlun 
experlnents. At lower temperatures this comparison Is obscured 
when the olivine becomes largely overgrown with pyroxene and/or 
groundmass . 

A restarkable feature of the ability of olivine to maintain 
surface equilibrium with the residual liquid. Is the constancy of 
the llquldus olivine composition. Irrespective of the cooling 
rate. The most magnesian olivines observed are a /out Fo^^. 
(Soxietlaes slightly less magnesian olivines are found, presumably 
because of failure to section the olivine core.) This coaposltlon 
Is observed Irrespective of the cooling rate and can be seen In 
Figure 7 to be the equilibrium llquldus olivine. Obviously the 
Roeder-EnsHe fomuiatlon of Fe/!4g distribution between 
olivine and liquid Is quite general and Independent of cooling 
rate over this range of cooling rates. 

This tendency of olivine to maintain surface Fe/Mg equilib- 
rium with the liquid even at very high cooling rates Is remarkable 
and was not expected. A growing olivine crystal must preferentially 
Incorporate Mg, and to a lesser extent Fe, while rejecting the 
other constituents of the liquid. Diffusive transport of material 
must then occur near the surface of the growing crystal, and at 



fast coollns rates eoaposltlon gradients in the liquid sight be 

expected to preclude equillbriua partitioning of Pe and Mg. A 

clue to the oeehanlsa for achieving such equillbriua may lie in 

# 

the skeletal nature of the blivines foraed at high cooling rates. 
Under these conditions any protrusions out froa the olivine growth 
surface aay penetrate the depleted liquid tone and thus saaple 
freah liquid without the need for long-range diffuaion. Such 
protrusions will therefore propagate rapidly yielding a skeletal 
crystal with near-equilibrlua Pe/Hg ratio. 

COCLINO hlSTORT 

We have seen that texture and crystal aorphology change as 
a function of a saaple *s cooling rate. Likewise pyroxene cheaistry 
shows soae variation with cooling rate. We have already noted . 
the qualitative nature of these variations and the crudeness of 
estimates of cooling rate based solely on such criteria. Clearly 
these criteria arc not sufficient for constructing a quantitative 
cooling history. They aay be used to check for consistency but 
they cannot point unaobiguously to soae particular cooling history. 
To interpret a cooling history quantitatively requires the ability 
to specify cooling rate for as oany stages in the crystallization 
history as possible. 

Since we do not know the possible effects of large super- 
heating, the earliest stage of the cooling history we can hope 
to reconstruct is the crystallization of olivine, the liquldus 
phase. Crystal aorphology is a rough guide but not really 


quantitative. The chenlstry of olivine appears to be insensitive 
to cooling rate and so ia no help in this undertaking. 

We have noted that the coaposition of the residual liquid 
at any teciperature appears to be roughly Independent of the 
cooling rate by which that temperature is reached. For this 
to be true, the total aaount of olivine crystallized prior to 
reaching a given temperature aust also be independent of cooling 
rate. This appears to be true in our experlnental runs; however, 
at slow cooling rates olivine occurs as relatively sparse etjuant 
crystals of relatively large volume, while at fast cooling rates 
olivine occurs as abundant thin plates of sxtall individual voIuxm. 
These observations suggest that the processes of crystal nuclcation 
and crystal growth may play reciprocal roles as the cooling rate 
is varied. In order to evaluate this possibility we counted the 
number of olivine crystals per unit area (a steasure of nucleatlon 
density) in our experimental products. The results are shown as 
a function of cooling rate in Figure 9> In spite of some scatter, 
a regular pattern is obvious; slow cooling produces a snail 
number of relatively large equant crystals, faster cooling producea 
many skeletal plates of small volume. Increasing the cooling rate 
for 30*/hour to A/1000*/hour does not result in further increase 
in nucleatlon density. Presumably, still faster cooling rates 
would result in reduced numbers of crystals and more quenched 
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liquid. Cooling rates faster than 10‘^*C/hour are sufficiently 
rapid to produce glasses with no olivine crystals. The nuclcation 
density, therefore, is at a maximum at cooling rates of a few 
hundred degrees per hour. Each block on Figure 9 includes runs 


^u«nch«d free Mny different tenpereturee at the Indicated cooling 
rate; the nuDber of crystala doea not vary ayateaatlcally with 
teeperature froa which the run was quenched, Isplylng that aost 
of the nuclei are formed early In the cooling process and that the 
observed scatter Is largely due to other factors. One such factor 
Is the snail sample else; another may be non-random distribution 
of olivine crystals with respect to the plane of the section. 

Some minor effect of continued nucleatlon of olivine with falling 
temperature Is also possible. In view of the large range of 
nucleatlon densities observed (1 crystal/surface to 200 crystals/ 
surface), the scatter does not preclude use of this paraaeter as 
an estlxator of cooling rates. 

The nunber of olivine crystals per unit area In thln-sectlon 
12002,9 Is also shown In Figure 9; a cooling rate of l*/hour, or 
a little less. Is Indicated for the beginning of olivine crystalli- 
sation In 12002. This estimate Is compatible with our earlier 
estimate baaed on olivine crystal morphology. 

This approach cannot be applied directly to estimate the 
cooling rate at the time of pyroxene crystallization. Since the 
beginning of pyroxene crystallization, unlike olivine crystalli- 
zation, occurs at progressively lower temperatures as cooling 
rate Is increased, the liquid from which pyroxene begins to 
crystallize varies In composition at different cooling rates. If 
12002, as we believe, did not cool at a uniform rate, estimates 
based on nucleatlon densities In our experimental runs might be 
misleading. As we have noted, the low Tl-^Al content of the 12002 


pyroxene phenocrysts compared to our experimental products 
(Figure 6) suggests that during this stage of Its cooling history, 
12002 was cooling at a rate slower than any we observed experi- 
mentally. A similar conclusion Is Implied by the large size (up 
to 3 am) of the natural pyroxene phenocrysts compared to those 
grown In our experiments (less than 0.5 mm). 

This conclusion can be made roughly quantitative by observa- 
tion of one aspect of pyroxene crystal morphology. The prominent 
"swallowtail” growth Instabilities seen on the pyroxene pheno- 
crysts In 12002 and In our experiments appear to have a charac- 
teristic linear spacing (Dj^X* In our experliaents this spacing 
Is a function of cooling rate. Figure 10 Is a plot of In(Dj^) 
versus cooling rate. The largest spacing shown by a symmetrical 
"swallowtail”. In each section Is plotted. The scatter Is sttrlbuted 
to slight variations In position and orientation of the crystals 
relative to the polished section. The range of "swallowtail" 
dimensions seen In 12002,9 Is also shown. Extrapolation of the 
experimental values suggests that the natural phenocrysts grew 
at a rate of 0.1 or 0.2*/hour, consistent with the qualitative 
estlBiates and several times slower than the cooling rate Inferred 
for the beginning of olivine crystallization. 

The crystallization of 12002 ends with groundmass growth 
and this event Is also dependent upon the previous cooling 
history. Again nucleatlon density estimates from our experiments 
at linear cooling rates cannot be rigorously applied since It 
appears that the rock did not cool at a constant rate. Furthermon 
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tht groundMss It not • tingle mineral but teveral mlneralt which 
grow penecontemporaneoutly . An additional complication la that 
nucleatlon density (or Itt equivalent here, grain tlxe) la not 
constant even at a tingle cooling rate. For example. Figure AC 
shows two regions of groundaass with different degrees of coarse- 
ness. Likewise there Is quite a variation In the grain site within 
the 12002 groundaass; a few patches have a coarseness equivalent 
to the coarsest seen In our experlaents even though the bulk of 
the 12002 groundaass Is much coarser (see Figure 11). While a 
quantitative estimate of the cooling rate of the 12002 groundaass 
Is not possible. It Is quite easy to tell that the average ground- 
aass of 12002 Is very auch coarser than any of our experimental 
products (as much as a factor of 10 or aore In crystal size). 

Since we have observed experlmqntally a coarsening of grain size 
In the groundsiass with decreasing cooling rate, the obvious con- 
clusion Is that the 12002 groundmass recoz*ds a slower cooling 
rate than any of our experlaents. 

12002 started crystallization with olivine, at an experi- 
mentally observed rate and concluded crystallization with ground- 
cass at a rate slower than any observed. Hence this porphyry 
solidified at a cooling rate that decreased with time. This 
contrasts with the cooling history we had originally Inferred 
froa the crystallization history recorded by porphyrltlc texture 
and Tl/Al variation In pyroxenes (Grove and others, 1973)* These 
considerations are Illustrated In Figure 12. The time scale on 
the dashed curve is arbitrary; however, the solid curve of the 


present Interpretation makes as much use of t.ie quantitative 
deductions as possible. The slope (cooling rate) at the start 
is the result of the estimate from olivine nucleatlon density. 

The control on the slope Is less firm at the temperatures corres- 
ponding to pyroxene and groundmass entry, but has been drawn to 
be roughly commensurate with possible estimates based on pyroxene 
Instability spacing and on the coarseness of the groundmass. 

Since the lower slopes are not well known the estimates of the 
time for the completion of solidification are somewhat elastic, 
although much of the temperature drop leading to solidification 
must take place In the first 20 days. The estimate for solidifi- 
cation time shown here Is on the order of a month or two, but 
It should be noted that It really Is pointless to try to specify 
this number more precisely since 12002 obviously never completely 
crystallized (there Is mesostasls glass) and the termination of 
the solidification was not a discrete event. What might be useful 
to attempt, though, would be refinements of the estimate of the 
cooling rate In the temperature Interval below the appearance of 
the groundmass. This might be attempted through cation ordering 
studies although It Is doubtful that sufficient calibration 
exists to substantially improve the quantitative resolution of 
the present conclusions. These conclusions are that 12002 began 
cooling at r* l*C/hour and that this cooling rate subsequently 
decreased. 

The form of the cooling history curve shown In Figure 12 Is 
npt consistent with the alternate Interpretation shown In which 


phcnocrysts grow during slow cooling (perhaps at depth) and 
grourdiiass grows during rapid quenching (perhaps after eruption). 
The deduced cooling history curve corresponds to that expected 
for slsple conductive cooling of magma solidifying under the 
influence of heat loss through neighboring margins of a ma^na 
body. The general cooling regime observed in solidifying terrest- 
rial Hawaiian lava lakes (Peck and others, 1966; Richter and 
Moore, 1966), is also for a decreasing cooling rate as indicated 
by the square root of time dependence of the thickening of the 
upper crust. However, in cooling through the melting interval 
in these lava lakes the cooling rate may Initially increase before 
decreasing (T.L. Wright, data on Makaopuhi lava lake, in press). 
Such effects are probably related to specific boundary conditions 
such as crust formation which do not seem to have influenced the 
cooling history of 12002. 

Conclusions 

We have shown here that finite cooling rates can radically 
alter the temperature and order of appearance of the crystalliza- 
tion phases from a silicate liquid. Using the results of 
experiments we have shown that the porphyritic texture of 12002 
could have formed in a single continuous cooling event. Nucleation 
density, swallowtail spaclngs, and groundaass grain size show 
that crystallization occurredat a decreasing cooling rate in a 
period of one or two months. 
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TABLE 1 

BUB rMOOUCTS fOt 12002 EOUILZBRATCO ZB NIQR 
PUBITT I BOB IB EVACUATCO SILICA TUBES 



#180 

#179 

#188 

•184 

•178 

#179 

#177 

#178 

#168 


LIQUID 

LIQUID 

LIQUID 

LIQUID 


UOUID 

LIQUID 

LIQUID 

LIQUID 

SiO} 

41. a 

41.7 

41.1 

44.0 

46.0 

46.0 

.49.8 

49.8 

44.9 

TlO] 

2. BO 

2.88 

2.41 

1.11 

1.91 

1.79 

1.78 

4.10 

4.74 

A1,0, 

7.B5 

8.22 

8.78 • 

9.81 

10.2 

11.1 

11.4 

11.1 

10.8 

cr,o. 

O.BS 

0.81 



0.97 

0.49 

0.18 

0.18 

0.18 

two 

22.0 

21.8 

22.8 

21.1 

19.9 

18.8 

19.1 

19.8 

19.0 

N90 

11. S 

12.9 

9.91 

8.89 

7.71 

7.00 

8.41 

8.28 

8.01 

HflO 

0.2« 

0.19 



0.12 

0.12 

0.11 

0.12 

0.18 

CaO 

•.IS 

8.19 

9.11 

9.99 

10.4 

11.2 

11.2 

11.6 

11.9 

Ba^O 

0.22 

0.07 



0.18 

0.04 

0.14 

0.14 

0.21 


•T.19 

90.74 

98.19 

98.90 

98.41 

90.70 

98.90 

99.18 

98.72 



Ota* 

OLZ- 

ota- 

OLI- 

oia- 

oca- 

(xa- 

OLI- 


VIBE 

:OBE 


VliS 


vns 

VIBE 

VI BE 

SIO, 

18.0 

19.9 

19.9 

18.8 

18.4 

19.9 

18.1 

19.8 

TlO] 

0.14 

O.ll 

0.09 

0.17 

0.19 

0.19 

0.20 

0.22 

CriO, 

0.78 



0.68 


0.74 

0.98 

0.49 

rao 

21.8 

29.8 

28.4 

27.8 

29.2 • 

M.8 

11.4 

12.0 

"90 

M.9 

19.9 

14.0 

11.8 

12.2 

11.0 

10.8 

10.1 

MAO 

0.29 



0.18 


0.40 

0.18 

0.40 

CaO 

0.14 

0.40 

0.19 

0.41 

0.41 

0.42 

0.49 

0.47 


99.89 

97.71 

98. 78 

99.80 

98.18 

99.49 

99.89 

99.68 


♦PTBO- 4PYBO- «BYBO- 9PYBO- 
EEKCS XEMES XEBES XEBES 


sto. 






91.7 

91.1 

92.9 

91.7 

TiO] 






0.98 

0.62 

0.70 

1.09 

Al2©3 






1.19 

1.29 

1.99 

2.21 

Cr,Oi 






0.79 

0.71 

0.78 

0.81 

Pao 






17.1 

17.8 

17.7 

18.1 

"90 






22.9 
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21.6 

20.7 

"no 







0. 18 

0.19 

0.19 

CaO 






1.80 

4.29 

4.92 

9. 99 

BajO 






0.01 

0.02 

0.09 

0.02 







98.01 

97.79 

100.19 

100.91 

Othar 

Phaaaa 



Crap 

Crap 

Crap 

Crap 

Crap 

Crap 

Prasar 

tt 







Plag 

Plag 

Bra Bua 









Dura* 

27 

28.1 

19.1 

18.7 

47 

72.7 

49.8 

98.7 

79 

tlon 










T*C 

1190 

1128 

1288 

1212 

1201 

1178 

1184 

1198 

1190 

Botaa 



Nlaor 

Minor 









Pa «aln 

Pa 9ala 







B**fl #171 #4 1104 C for 1B3«4 boors hod sassobls^o oIIvIao, pyroxono, plso 
•p. lloonlta. 
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FIGURE CAPTIONS 

Fieuf 1 

Fhas* relations for 1?002 bulk eosiposltlon deternlned by 
syntheala experiaenta on dcvltrlfied glaaa atartlng material. 
Pleura 2 

Reaulta of controlled cooling experlmenta as a function of 
the cooling rate and temperature of quenching. Dotted lines 
Indicate cooling trajectories froai l350t$*C, where all experiments 
start, to the quenching temperatux^. All experiownts below a 
particular curve, such as "olivine in", have olivine present in 
the product. These phase appearance curves nay be compared to the 
equilibrium crystallixation of 12002 shown in the upper right 
portion of the figure. 

Figure 3 

Products of rapid cr-oling. Both photomicxx>graphs to same 
scale. Bar is 100 pm. Upper photo is a transmitted light view 
of thin section 12009 seen with crossed nicols. Olivine forms 
long skeletal blades with elaborate cross-ornamentation. In the 
center of the picture is a cluster of dendritic pyroxene. Black 
areas are glass and some subaicroscoplc devitrification products. 
The lower photograph is a reflected light view of an experiment on 
12002 cocpositlon cooled at 'v500*C/hour. Note the sisiilarities to 
the upper photo: barred olivine blades, a cluster of dendritic 
pyroxene in the left-center, and glass matrix. Bubbles on bottom 
and top- left are in the epoxy mounting medium. 

Figure A 

Experimental textures seen in reflected light photomicro- 
graphs. Scale the same on each photo; bar is 100 pa. Cooling 


rate decreases from A(57*C/hour) , B(ll*C/hour) , C(3*C/hour), to 
D(2*C/hour). ^ Skeletal olivine with spinel specks around edges. 
Skeletal high-Ca,Al pyroxenes, llaenite fibers resolvable in glass 
groundmass. B^ Principle change from A is that pyroxenes are 
more euhedral and groundmass is microcrystalline plag^px (♦tsilica). 
Olivine is overgrown by pyroxene. Ilmenite is easily observable. 

C. Pyroxenes and olivines similar to B. Groundmass is coarser 
and in right center individual crystals are resolvable. Note that 
the groundmass is not this coarse in all portions of the surface. 

D. Olivine and pyroxene somewhat coarser than C but groundmass 
is much coarser than in C. 

Figure 5 

Spinel relations — a study in heterogeneous nucleation. 

Reflected light photomicrographs of experimental and natural spinels 
Scale bar is 30 pm in each photo. ^ Minute spinel specks orna- 
menting the interior trellis of a skeletal olivine. Glass forms 
the background. Grown at 1200*C/hour. Dendritic spinel cross 

growth on the margin of an olivine phenocryst (lower portion of 
photo). Skeletal high-Ca,Al pyroxenes and glass fill upper portion 
of photo. Grown at 57*C/hour. Mantled spinel phenocryst. Cr- 

spinel interior and ulvospinel rim have grown around an iron globule 
A bubble appears to have been trapped. Coarse groundmass with 
ilmenite laths is seen to the left between olivine and pyroxene 
phenocrysts. Crown at 3*C/hour. D^ Natural mantled spinels in 
12002,9. Note similarities to C of the association with iron. 

Figure 6 

Compositional distribution of pyroxenes observed in 12002, 
controlled cooling experiments, and equilibrium crystallisation of 



12002. Data for 12002 ara froa Crova and othara (1973). Exparl- 
aantal pyroxanaa ara all producad in Iron capsulaa In aaalad 
silica tubas. 

7 

Fa/Fa^M£ of 12002 axparlawntal ollvlnas aqd liquids Is shown 
as a function of taaparatura for axparlaants In Iron capsulaa In 
saalad silica tubas. Solid clrclas and trlan^las ara for ollvlnas 
and liquids In aqulllbrlua axparljsants. Bars Indlcata the range 
In olivine cocpos It Ion observed In controlled cooling runs. Bars 
ara shown at the quenching tanparatura of the axparlaent, but the 
range to higher Po contents Indicates that the ollvlnas record 
the high tarparatura parts of the cooling exparlnant. Open triangles 
ara liquid coaposltlons in controlled cooling runs. Compositional 
gradients ara observed In the liquid at high cooling rates and at 
low tanparaturas. For Inhomogeneous liquids, the open triangles 
record the liquid composition adjacent to the olivine rims of 
highest Fe/Fe-^Mg. Kaxlmum Fe/Fe-t-Mg In ollvlnas and liquids In 
the controlled cooling runs approximate the equilibrium values as 
a function of tenperature above the equilibrium solidus. Deviations 
do not correlate with cooling rate. Minimum Fe/Fe-fNg values 
observed In the ollvlnas are shown on a histogram In the lower 
left. the equilibrium llquldus olivine. Is most often 

observed, Is observed at all cooling rates and the deviations 

are not correlated with cooling rate. Low Fe/Fe^^Mg deviations are 
thought to be the result of iron capsule impurities. High Fe/Fe-fNg 
deviations nay result from oxidation of the capsule during prepara- 
tion or from failure to section or observe the Bw>st forsterltle 
olivine core. 


Figure 6 

Ollvlne-llquld partitioning of Fe and Kg. Data from Figure 7 
for olivines (maximum Pe/Fe*Mg observed) and liquids are compared 
to curves for distribution coefficients (X^ of Boeder and Emslle, 
1970). Equilibrium experiments show very little scatter about 
the curve Kp * 0.33, whereas controlled cooling experiments show 
some scatter. This scatter Is not great and does not cox*relate 
with cooling rate. 

Figure 9 

Nucleatlon density of olivine vs cooling rate In experimental 
charges. The natural logarithm of the number of olivine crystals 
observed per mm^ of section 'surface Is used as a measure of nuclea- 
tlon density. This parameter has been measured In thin section 
12002,9 and yields an estimate slightly less than l*C/hour for the 
cooling rate* at the beginning of 12002 crystallisation. 

Figure 10 

The natural logarithm of the spacing of growth Instabilities 

(D,..) In pyroxene Is plotted versus the natural logarithm of the 
13 

cooling rate at which the pyroxenes grew. D^^ Is the half width 
of "swallow tails." The range in this paraiseter for 12002,9 
pyroxene phenocrysts is shown and an extrapolation Indicates that 
they may have grown at a cooling rate of 0.1-0. 2* C/hour. 

Figure 11 

Groundmass relations shown In section 12002,9. Scale bar Is 
100 urn In all pictures. A and A* (reflected and transmitted light) 
show a fine-grained fascicle of pyroxene and plagloclase. The 
pyroxene Is the more highly reflective. The grain slse In this 
fascicle Is comparable to the coarsest seen In our experiments 
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(see Figure %0) but Is styplesl of the major portion of the 12002 
grounOasss which Is coarser. Note coarser texture surrounding 
this fascicle. B and B' show a coarser fascicle than A and A', 
but again notice a small, fine-grained fascicle In the upper right. 
Also seen are a variety of llsienlte habits. C and 0 Illustrate 
the difficulties In determining the place of llsMnlte In the 
crystallisation sequence on the basis of textural relations. This 
is also Illustrated In B. Ilmenlte precedes groundsiass crystalli- 
sation In linear cooling experlswnts and may do so In 12002 also 
but the textural evidence Is asd>lguous. 

Figure 12 

Susaary of present results compared to our previous Interpre- 
tation (Grove and others, 1973) • Present conclusions are based on 
comparison of 12002 with controlled cooling experiments. Our 
previous conclusions were based on a comparison of 12002 with 
equilibrium experiments. (Time scale on this dashed curve Is 
srbltrary . ) 
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